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Abstract 

We extend previous computations of Calabi-Yau metrics on projective 
hypersurfaces to free quotients, complete intersections, and free quotients 
of complete intersections. In particular, we construct these metrics on 
generic quintics, four-generation quotients of the quintic, Schoen Calabi- 
Yau complete intersections and the quotient of a Schoen manifold with 
Z 3 x Z 3 fundamental group that was previously used to construct a het- 
erotic standard model. Various numerical investigations into the depen- 
dence of Donaldson's algorithm on the integration scheme, as well as on 
the Kahler and complex structure moduli, are also performed. 
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1 Introduction 



A central problem of string theory is to find compactifications which can reproduce 
real world physics, in particular the Standard Model. The first and still one of the 
best motivated ways to achieve this are heterotic string compactifications on Calabi- 
Yau manifolds pQ. In particular, the so-called "non-standard embedding" of E$ x E$ 
heterotic strings has been a very fruitful approach towards model building. 

For a variety of reasons, the most successful models of this type to date are based on 
non-simply connected Calabi-Yau threefolds. These manifolds admit discrete Wilson 
lines which, together with a non-flat vector bundle, play an important role in breaking 
the heterotic Eg gauge theory down to the Standard Model [21 El HI El El [7j . In the 
process, they project out many unwanted matter components. In particular, one can 
use this mechanism to solve the doublet-triplet splitting problem [H [9] . Finally, the 
non-simply connected threefolds have many fewer moduli as compared to their simply 
connected covering spaces [10]. In recent work [HI [T2J, [131 H3] , three generation models 
with a variety of desirable features were introduced. These are based on a certain 
quotient of the Schoen Calabi-Yau threefold, which yields a non-simply connected 
Calabi-Yau manifold with fundamental group Z3 x Z3. 

Ultimately, it would be desirable to compute all of the observable quantities of 
particle physics, in particular gauge and Yukawa couplings, from the microscopic 
physics of string theory [T5l [T6] . There are many issues which must be addressed to 
do this. Physical Yukawa couplings, for example, depend on both coefficients in the 
superpotential and the explicit form of the Kahler potential. In a very limited number 
of specific geometries [TTJ, [HI [191 I2H] , the former can be computed using sophisticated 
methods of algebraic geometry, topological string theory and the like. For the latter, 
we generally have only the qualitative statement that a coefficient is "expected to be 
of order one" . Doing better and, in particular, extending these calculations to non- 
standard embedding, multiply connected compactifications has been an outstanding 
problem [I]. 

In recent work [2U [22] , a plan has been outlined to analyze these problems numer- 
ically, at least in the classical limit. The essential point is that, today, there are good 
enough algorithms and fast enough computers to calculate Ricci-flat metrics and to 
solve the hermitian Yang-Mills equation for the gauge connection directly. Given this 
data, one can then find the correctly normalized zero modes of fields, determine the co- 
efficients in the superpotential and compute the explicit form of the Kahler potential. 
Some progress in this direction was made in [2H [221 (231 [2H [22] , and also (23, [23 [21] • 

In the present work, we take some significant steps forward in computing Calabi- 
Yau metrics. Some of the steps are technical and computational improvements, which 
we will discuss. But the primary new ingredient is the ability to solve for Ricci-flat 
metrics on non-simply connected Calabi-Yau manifolds. While in broad conceptual 
terms the procedure is similar to the simply connected case, in practice the problem of 
finding and working with a complete basis of holomorphic sections of a line bundle, as 
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used in Donaldson's method, is now quite intricate. To solve this, we systematically 
use the techniques of Invariant Theory [29] . 

We begin, in ISection 2| by extending Donaldson's algorithm to the computation 
of Calabi-Yau metrics for generic quintic threefolds. This formalism is first applied 
to the simple Fermat quintic, reproducing and extending the results of [23]. We 
then numerically calculate the Calabi-Yau metrics, and test their Ricci-flatness, for 
a number of random points in the complex structure moduli space. All of these 
manifolds are, of course, simply connected. We then proceed to non-simply connected 
manifolds or, equivalently, to covering spaces that admit fixed point free group actions. 
In ISection 3| we outline the general idea and review some of the Invariant Theory, 
in particular the Poincare series, Molien formula and the Hironaka decomposition, 
that we will use. This formalism will then be applied in ISection 4l to the subset of 
quintics that admit a Z5 x Z5 fixed point free group action. Using the Molien formula 
and the Hironaka decomposition, we determine the Z5 x Z5 invariant sections on such 
quintics and, hence, on the Z 5 x Z 5 multiply connected quotient space. Using these, 
we can extend Donaldson's algorithm and compute the Calabi-Yau metric, and test 
its Ricci-flatness, on the quotient. As a by-product of this process, we note that there 
are now two ways to compute a Calabi-Yau metric on the covering space; first, as a 
point in the quintic moduli space employing the methods of ISection 21 and second, 
by using Donaldson's algorithm for invariant sections only. These two methods are 
compared in ISection 4l Note that only the second approach descends to the Z 5 x Z 5 
quotient threefolds. 

In ISection 5l we describe Schoen threefolds. We show how to compute Calabi-Yau 
metrics on these simply connected complete intersection manifolds and, as always, test 
their Ricci-flatness. Schoen manifolds which admit a fixed point free Z3 x Z3 group 
action are then discussed in ISection 61 Proceeding as in ISection 4| the Molien formula 
and the Hironaka decomposition are generalized to complete intersections. These are 
then used to find all Z 3 x Z 3 invariant sections. These descend to the Z 3 x Z 3 quotient, 
and are used to compute the Calabi-Yau metric on this multiply connected threefold. 

In addition, we explicitly check algebraic independence of primary invariants, 
which are defined in ISection 3l for quintics in Appendix A Finally some portions 



of the code used in this paper are presented in Appendix B 



2 The Quintic 

2.1 Parametrizing Metrics 

Quintics are Calabi-Yau threefolds Q C P 4 . As usual, the five homogeneous coordi- 
nates [z : Z\ : z 2 ■ z 3 : Z4] on P 4 are subject to the identification 

[zq : Z\ : z 2 : z 3 : z 4 ] = [Xz : Xzi : Xz 2 : Az 3 : Az 4 ] VA G C — {0}. (1) 
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In general, a hypersurface in P 4 is Calabi-Yau if and only if it is the zero locus of a 
degree-5 homogeneous polynomials 



Q( Z ) — C (n ,n 1 ,n 2 ,n3,n4) z 0° z ± 1 z 2 2 Z 3 3 ■ (2) 

"0+™i+«2+n.3+"4=5 

Note that, abusing notation, we denote both the threefold and its defining poly- 
nomial by Q. There are ( 5+ 4~ 1 ) = 126 degree-5 monomials, leading to 126 coef- 
ficients C( no>nijn2in3in4 ) G C. These can be reduced by redefining the ^-coordinates 
under GL(5,C). Hence, the number of complex structure moduli of a generic Q is 
126 — 25 = 101. A particularly simple point in this moduli space is the so-called 
Fermat quintic Qf, defined as the zero-locus of 

Q F (z) = z\ + z\ + 4 + 4 + 4- ( 3 ) 

We will return to the Fermat quintic later in this section. 

In general, the metric on a real six- dimensional manifold is a symmetric two- 
index tensor, having 21 independent components. However, on a Calabi-Yau (more 
generally, a Kahler) manifold the metric has fewer independent components. First, 
in complex coordinates, the completely holomorphic and completely anti-holomorphic 
components vanish, 

gij(z,z)=0, gy(z,z) = 0. (4) 
Second, the mixed components are the derivatives of a single function 

9is(z, z ) = z ) = didjK{z, z). (5) 

The hermitian metric g^ suggests the following definition of a real (1, l)-form, the 
Kahler form 

% % — 

to = -g rj dzi A dz- 3 = -ddK(z, z). (6) 

The Kahler potential K (z, z) is locally a real function, but not globally; on the overlap 
of coordinate charts one has to patch it together by Kahler transformations 

K(z,z)~K{z,z) + f(z) + f{z). (7) 

The metric eq. (jSJ) is then globally defined. 

The 5 homogeneous coordinates on P 4 clearly come with a natural 577(5) action, 
so a naive ansatz for the Kahler potential would be invariant under this symmetry. 
However, the obvious S77(5)-invariant \zq\ 2 + ■ ■ ■ + 1 2r 4 1 2 would not transform correctly 

1 Hence the name quintic. 
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under the rescaling eq. (CD) with A = X(z). Therefore, one is led to the uniqu^l SU(5) 
invariant Kahler potential 

l 4 

K FS = - In yZiZj. (9) 

IT < * 



7T 

i=0 



One can slightly generalize this by inserting an arbitrary hermitian 5x5 matrix h a @, 

i 4 

K FS = -ln V h a(3 z a z . (10) 
7r ^rr" 

a,/3=0 

Any Kahler potential of this form is called a Fubini-Study Kahler potential (giving rise 
to a Fubini-Study metric). At this point the introduction of an arbitrary hermitian 
h a/3 does not yield anything really new, as one can always diagonalize it by coordinate 
changes. However, strictly speaking, different h a/3 are different Fubini-Study metrics. 

The above Kahler potential is defined on the whole P 4 and, hence, defines a metric 
on P . But this induces a metric on the hypersurface Q, whose Kahler potential is 
simply the restriction. Unfortunately, the restriction of the Fubini-Study metric to 
the quintic is far from Ricci-flat. Indeed, not a single Ricci-flat metric on any proper 
Calabi-Yau threefold is known. One of the reasons is that proper Calabi-Yau metrics 
have no continuous isometries, so it is inherently difficult to write one down analyt- 
ically. Recently, Donaldson presented an algorithm for numerically approximating 
Calabi-Yau metrics to any desired degree ^ZT\. To do this in the quintic context, take 
a "suitable" generalization, that is, one containing many more free parameters of the 
Fubini-Study metric derived from eq. ffTO]) on P 4 . Then restrict this ansatz to Q and 
numerically adjust the parameters so as to approach the Calabi-Yau metric. An ob- 
vious idea to implement this is to replace the degree- 1 monomials z a in eq. ( flOl) by 
higher degree-/? monomials, thus introducing many more coefficients in the process. 
However, note that the degree k is the Kahler class 

k G H 1 ' 1 ^ 4 , Z) ~ Z. (11) 

The reason for this is clear, for example, if we multiply Kps in eq. ([9]) by k. Then 



k 



4 1 4 

In^ZjZj = - In ^2 z h ' ' - z i k Zn ■ ■ ' z %- ( 12 ) 



FS - 

7T ' — ' 7T 

i=0 ii,...,ifc=0 



2 Uniquc up to an overall scale, of course. The scale is fixed by demanding that wfs is an integral 
class, lo G i? 2 (P 4 ,Z). To verify the integrality, observe that the volume integral over the curve 
[1 : t : : : 0] in P 4 is 

/ % -ddK FS ([l : t : : : 0]) = / -«9 t <% M 1 + ^ dt d * = L ( 8 ) 



5 



k 


1 


2 


3 


4 


5 


6 


7 


8 




5 


15 


35 


70 


126 


210 


330 


495 


N k 


5 


15 


35 


70 


125 


205 


315 


460 



Table 1: The number of homogeneous polynomials N k and the number of re- 
maining polynomials N k after imposing the hypersurface constraint, 
see eq. (fT6j) , 



Hence, if we want to keep the overall volume fixed, the correctly normalized general- 
ization of eq. (TIOl is 



I 4 

K(z, z) — — In h^'-M'&'-M z h --- z ik % • • • z % . (13) 



kir 

h,...,ik=0 
3l,-,Jk=° 



degree k degree k 



Note that the monomials z^ ■ ■ ■ z ik , where z'i, . . . , = 0, . . . , 4 and integer k > 0, are 
basis vectors for the space of polynomials C[zq, . . . , z k }. For fixed total degree k, they 
span the subspace C[z , . . . , z k ] k of dimension 

Some values of N k are given in lTablell In particular, the matrix of coefficients h now 
must be a hermitian Nk X Nk matrix. 

However, there is one remaining issue as soon as k > 5, namely, that the monomials 
will not necessarily remain independent when restricted to Q. In order to correctly 
parametrize the degrees of freedom on Q, we have to pick a basis for the quotient 

C[z ,...,z 4 ] k /(Q(z)) (15) 

for the degree-/? polynomials modulo the hypersurface equation. Let us denote this 
basis by s a , a — 0, . . . , — 1. It can be shown that for any quintic 

< k < 5 

N k = { " (16) 

>-^=r*- i )-(a k>5. 

Some values of the Nk are listed in ITableTl For any given quintic polynomial Q(z) 
and degree k, computing an explicit polynomial basis {s a } is straightforward. As an 
example, let us consider the the Fermat quintic defined by the vanishing of Qf{z), see 
eq. Q- In this basis for the quotient eq. f[T5l) can be found by eliminating from 

any polynomial in C[zq, . . . , z^ k all occurrences of z\ using z\ = — (z\ + z\ + z\ + 
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Using the basis s a for the quotient ring, one finally arrives at the following ansatz 

N k -1 

a,/3=0 

for the Kahler potential and, hence, the approximating metric. Note that they are 
formally defined on P 4 but restrict directly to Q, by construction. Obviously, this 
is not the only possible ansatz for the approximating metric, and the reason for this 
particular choice will only become clear later on. However, let us simply mention 
here that there is a rather simple iteration scheme [23l [2U [30] involving which 
will converge to the Ricci-flat metric in the limit k — > oo. Note that, in contrast 
to the Fubini-Study Kahler potential eq. (fill , the matrix h a/S in eq. ( fTTj) cannot be 
diagonalized by a GL(5, C) coordinate change on the ambient P 4 for k > 2. 

Let us note that there is a geometric interpretation of the homogeneous polyno- 
mials, which will be important later on. Due to the rescaling ambiguity eq. ([!]), the 
homogeneous polynomials are not functions on P 4 , but need to be interpreted as sec- 
tions of a line bundle. The line bundle for degree-/c polynomials is denoted P 4(/c) 
and, in particular, the homogeneous coordinates are sections of P 4(1). In general, 
the following are the same 

• Homogeneous polynomials of degree k in n variables. 

• Sections of the line bundle ¥ n-\{k). 

Moreover, the quotient of the homogeneous polynomials by the quintic, eq. ( fl5i) . is 
geometrically the restriction of the line bundle P 4(fc) to the quintic hypersurface. 
That is, start with the identification above, 

#°(P 4 ,(V(£;)) =C[z ,z 1 ,z 2 ,z 3 ,z 4 \ k . (18) 

After restricting the sections of P 4.(k) to Q, they satisfy the relation Q(z) = 0. Hence, 
the restriction is 

H°(Q,On(k)) =C[z Q ,z 1 ,z 2 ,z 3 ,z i ] k / (q(z)) 

')„ /k x (19) 

= C[Z0) Z\, Z2, Zzi Z^[k I ( QC[Zo, Z\, Z2, Zs, Z4]k-5 J ■ 

More technically, this whole discussion can be represented by the short exact sequence 
. H° ( P 4 , P 4 (k - 5)) X S- H° ( P 4 , p4 (jfe)) restrict > H° (Q, O^k)) 



z4k-5 XQ{Z) ' C [«0» • • • > Z i\k C N, ■ ■ ■ , Z A ] k / (Q(z)\ 



k 



(20) 
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2.2 Donaldson's Algorithm 

Once we have specified the form for the Kahler potential, our problem reduces to 
finding the "right" matrix h a ^. This leads us to the notion of T-map and balanced 
metrics, which we now introduce. First, note that eq. (fTTj) provides a way to define 
an inner product of two sections. While it makes sense to evaluate a function at a 
point, one cannot "evaluate" a section (a homogeneous polynomial) at a point since 
the result would only be valid up to an overall sea lfl However, after picking \\s 
see eq. (fTTl) . one can cancel the scaling ambiguity and define 



2 

h,ki 



{s ' SM= jw = zj^Fm Vp£Q (21) 

for arbitrary sections (degree-/c homogeneous polynomials) S, S' G H°(Q, Og(fc)). 
Note that the sq, . . . , Sjv fe -i are a basis for the space of sections, so there are always 
constants c a G C such that 

N k -1 

The point- wise hermitian form ( , ) is called a metric on the line bundle Og(k). Given 

this metric, we now integrate eq. ( T2TT) over the manifold Q to define a C-valued inner 
product of sections 



[S,S')= Nk [ (S, S') (p) dVolcY 

VoWfQ) Jo (23) 

dVolcY ■ 



icy(Q) 



Vol CY (Q) 7g Ea,P h^S a Sp 

Since ( , ) is again sesquilinear, it is uniquely determined by its value on the basis 
sections s a , that is, by the hermitian matrix 

H a - p = {s ai s p ). (24) 

In general, the matrices h a/3 and H a p are completely different. However, for special 
metrics, they might coincide: 

Definition 1. Suppose that 

h'*={H a$ y\ (25) 
Then the metric h on the line bundle Og{k) is called balanced. 

We note that, in the balanced case, one can find a new basis of sections {s a }a=o 1 
which simultaneously diagonalizes H a n = 5 a p and h a/3 = 5 a/3 . The interesting thing 
about balanced metrics is that they have special curvature properties, in particular 



3 In other words, at any given point one can only decide whether the section is zero or not zero. 
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Theorem 1 (Donaldson [30]). For each k > 1 the balanced metric h exists and is 
unique. As k oo, the sequence of metrics 



1 N k -l 
a,/3=0 

on Q converges to the unique Calabi- Yau metric for the given Kdhler class and complex 
structure. 

Hence, the problem of finding the Calabi- Yau metric boils down to finding the 
balanced metric for each k. Unfortunately, since H a § depends non-linearly on h a ^ 
one can not simply solve eq. (125]) defining the balanced condition. However, iterating 
eq. (125]) turns out to converge quickly. That is, let 



T(h) a g = H a s = — * I - a -§ _ dVolcv (27) 



Vol CY (Q) Jq Ztfh-rh 
be Donaldson's T-operator. Then 

Theorem 2 (Donaldson, [21]). For any initial metric h , the sequenc^ 

h n+1 = {T(h n )Y l (28) 
converges to the balanced metric as n — > oo . 

In practice, only very few (< 10) iterations are necessary to get very close to the 
fixed point. Henceforth, we will also refer to g~ in eq. (1261) . the approximating metric 
for fixed k, clS db balanced metric. 

2.3 Integrating over the Calabi- Yau threefold 

We still need to be able to integrate over the manifold in order to evaluate the T- 
operator. Luckily, we know the exact Calabi- Yau volume form, 

dVolcv = n a n, (29) 

since we can express the holomorphic volume form Q as a Griffiths residue. To do 
this, first note that the hypersurface Q C P has complex codimension one, so we can 
encircle any point in the transverse direction. The corresponding residue integral 

fl = I (30) 
/ Q{z) 



4 At this point, it is crucial to work with a basis of sections sq, . . . , spf k -i- For if there were a 
linear relation between them then the matrix T(h) would be singular. 
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is a nowhere vanishing holomorphic (3, 0)-form and, hence, must be the holomorphic 
volume form Q. As an example, consider the Fermat quintic defined by eq. ([3]). In 
a patch where we can use the homogeneous rescaling to set z = 1 and where z 2 , Z3, 
and Z4 are good local coordinates, 

/" dz 1 A ■ ■ ■ A dz4 _ dz 2 A dz 3 A dz 4 
~ J l + z% + z% + 4 + zl~ 5zf ' ^ ' 

To apply, for example, Simpson's rule to numerically integrate over the Calabi-Yau 
threefold one would need local coordinate charts. However, there is one integration 
scheme that avoids having to go into these details: approximate the integral by N. p 
random points {pi}, 

1 Np r 

TfE^)^ //dVol. (32) 

p i=l J 

Of course, we have to define which "random" distribution the points lie on, which 
in turn determines the integration measure dVol. In practice, we will only be able 
to generate points with the wrong random distribution, leading to some auxiliary 
distribution dA. However, one can trivially account for this by weighting the points 
with Wi = £7 A Q,/ dA, 



W r t = W p t ^jjVgu-ff dVo lcv . (33, 

Note that taking / = 1 implies 

1 Np 

^. = Vol CY (Q). (34) 



^ 1=1 



Points from Patches 

We start out with what is probably the most straightforward way to pick random 
points. This method only works on the Fermat quintic, to which we now restrict. Let 
us split P 4 into 5 • 4 = 20 closed sets 



- I [zo '■ Z\ : z 2 : Z3 : 24] 

\z e \ = max(|z |, • • • , N), \z m \ = max(|^ |, • • • , M, • • • , l^l)}- (35) 

In other words, Zg has the largest absolute value and z m has the second-largest absolute 
value. They intersect in real codimension-1 boundaries where the absolute values are 
the same and induce the decomposition 

Qf = [J Qf/tu (36) 
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with QF,£m = Qf H Ue m . Since permuting coordinates is a symmetry of the Fermat 
quintic, it suffices to consider Qf,oi- We define "random" points by 

• Pick x, y, z G C<i on the complex unit disk with the standard "flat" distribution. 

• Test whether 

\x\, \y\, \z\ < |l + x 5 + y 5 + z 5 \" < 1. (37) 

If this is not satisfied, start over and pick new x, y, and z. Eventually, the above 
inequality will be satisfied. 

• The "random" point is now 

1 : -(1 +x 5 +y 5 + z 5 ) J : x : y : z G Q Ffll , (38) 

where one chooses a uniformly random phase for the fifth root of unity. 

By construction, the auxiliary measure is then independent of the position (x, y, z) G 
Q F ,im- Hence, 

dA = — d 2 x AdflyA d 2 z. (39) 
20 y J 

Points From Intersecting Lines With The Quintic 

The previous definition only works on the Fermat quintic, but not on arbitrary quin- 
tics. A much better algorithm [23J is to pick random lines 



L-P'cP 4 . (40) 

Any line L determines 5 points by the intersection L H Q = {5pt.} whose coordinates 
can be found by solving a quintic polynomial (in one variable) numerically. Explicitly, 
a line can be defined by two distinct points 

P = [Po : Pi ■ P2 ■ P3 : Pi], q = [qo ■ qi ■ <?2 : <?3 : Qa] e P 4 (41) 

as 

L : C U {oo} -> P 4 , [p + tq : Pl + tq x : p 2 + tq 2 : p 3 + tq 3 : p 4 + (42) 

The 5 intersection points L D Q are then given by the 5 solutions of 

Q ° L{t) = Q(p + tq ,p! + tq 1 ,p 2 + tq 2 ,p 3 + tq 3 ,p 4 + tq A ) = 0. (43) 

Clearly, the auxiliary measure will depend on how we pick "random" lines. The 
easiest way is to choose lines uniformly distributed with respect to the SU{5) action on 
P 4 . Note that a line L is Poincare dual to a (3, 3)-current, that is, a (3, 3)-form whose 
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coefficients are delta-functions supported on the line L. For the expected distribution 
of lines, we then average over all "random" configurations of lines. Because of this 
averaging procedure, the Poincare duajfl of the expected distribution of lines (L) is a 
smooth (3, 3) form. Since there is (up to scale) only one 5'[/(5)-invariant (3, 3) form 
on P 4 , the expected distribution of lines must be 

(L) ~ c4 s , (44) 

where cups is the Kahler form defined by the unique S'?7(5)-invariant Fubini-Study 
Kahler potential eq. ([9]). Restricting both sides to an embedded quintic i : Q P 4 , 
we obtain the auxiliary measure as the expected distribution of the intersection points, 

dA= (QDL) ~r(c4 s ). (45) 

As a final remark, note that the the symmetry of the ambient space is, in general, not 
enough to unambiguously determine the auxiliary measure. It is, as we just saw, suf- 
ficient for any quintic hypersurface Q. However, for more complicated threefolds one 
needs a more general theory. We will have to come back to this point in lSubsection 5.41 



2.4 Results 

Following the algorithm laid out in this section, we can now compute the successive 
approximations to the Calabi-Yau metric on Q. In order to test the result, we need 
some kind of measure for how close the approximate metric is to the Calabi-Yau 
metric. Douglas et al. [23] proposed the following: First, remember that the Kahler 
form uj eq. ([6]) is the Calabi-Yau Kahler form if and only if its associated volume form 
u 3 is proportional to the Calabi-Yau volume form eq. fl30|) . That is, 

u 3 (p) = (const.) x (n(p) A n(p) ) Vp G Q (46) 

(the Monge-Ampere equation) with a non-vanishing proportionality constant inde- 
pendent of p E Cjfl Let us define 

Vol K (Q) = [ to 3 (47a) 

J Q 

and recall that 

VoIcy (Q) = [ ttAtt. (47b) 
J Q 

5 By the usual abuse of notation, we will not distinguish Poincare dual quantities in the following. 
6 And varying over the moduli space. However, this will not concern us here. 
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Points chosen from coordinate patches: a k 
Fit for k > 3: a { k l) = 3.23/ k 2 _- 4.55/fc 3 

Points from Q f] L: 
Fit for k > 3: af ] = 3A0/k 2 - 4.96 /k 3 



k = k = l k = 2 k = 3 k = A k = 5 k = 6 k = 7 k- 

d (k) 



Figure 1: The error measure c/% for the metric on the Fermat quintic, com- 
puted with the two different point generation algorithms described in 
\Subsection 2.31 In each case we iterated the T-operator 10 times, 
numerically integrating over N p = 200,000 points. Then we eval- 
uated <7fc using 10,000 different test points. The error bars are the 
numerical errors in the <7fc integral. 
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The ratio of these two constants determines the proportionality factor in eq. (T46 
This equation can now be rewritten 



Vol K (Q) Volov (Q) 

Note that one often demands that the two constants, eqns. (j47a|) and (14 7b |) . are unity 
by rescaling u and Q respectively. However, this would be cumbersome later on and 
we will not impose this normalization. Then the integral 



VoIcy (Q) Jq 



Vol K (Q) 



QAQ/ VoIcy (Q) 



dVolcv (49) 



vanishes if and only if u is the Calabi-Yau Kahler form. In practice, Donaldson's 
algorithm determines successive approximations to the Calabi-Yau metric. Since we 
know the exact Calabi-Yau volume form Q A Q, only u is approximate and depends 
on the degree k. We define to be the above integral evaluated with this degree-fc 
approximation to the Calabi-Yau Kahler form. 

Let us quickly summarize the steps necessary to compute the metric. To do that, 
one has to 

1. Choose a degree k at which to compute the balanced metric which will approx- 
imate the Calabi-Yau metric. 

2. Choose the number N p of points, and generate this many points {pi}f=i on Q. 
Although k and N p can be chosen independently, we will argue below that N p 
should be sufficiently larger than N% for accuracy. 

3. For each point p i: compute its weight Wi = dA(pi)/(Q A 0). 

4. Calculate a basis {s a }a=o 1 ^ or ^ ne quotient eq. fTT5|) at degree k. 

5. At each point calculate the (complex) numbers {saipi)}^^ 1 an d; hence, the 
integrand of the T-operator. 

6. Choose an initial invertible, hermitian matrix for h lS . Now perform the numer- 
ical integration 



nhu= f fc y (50) 



7. Set the new h a ' p to be h a ~ p = (T^)' 
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8. Return to litem 61 and repeat until h a ^ converges close to its fixed point. In 
practice, this procedure is insensitive to the initial choice of h a/S and fewer than 
10 iterations suffice. 

Having determined the balanced h a @, we can evaluate the metric using eq. (1261) 
and, hence, the Kahler form uj{p) at each point p, see eq. ((6]). Now form u> 3 (p). This 
lets us compute by the following steps: 



1. The o"fe integral requires much less accuracy, so one may pick a smaller number 
N p of poh 

2. Compute 



N p of points {pi}f= v 



Vol CY = ^J>, Vol K = i, £ " ^, (51) 

which numerically approximate Q A Q and Jq cj 3 , respectively. 
3. The numerical integral approximating Ok is 

^) 3 /Vol K 



iV P Vol CY 



E 



i=l 



Q{ Pi ) Afi(pi)/VolcY 



w». (52) 



As a first application, we apply this procedure to compute the Calabi-Yau metric 
for the simple Fermat quintic Qf defined by eq. (j3J). In this case, there are two point 
selection algorithms, both given in ISubsection 2.31 We do the calculation for each 



and show the results in Figure 1 One can immediately see that both point selection 
strategies give the same result, as they should. In fact, there is a theoretical prediction 
for how fast cr^ converges to 0, see [221 [23 ED] • Expanding in |, the error goes to zero 
at least as fast as 

°* = F + F + ""' 5iGR - (53) 

In particular, the coefficient of \ is proportional to the scalar curvature and vanishes 

& I S3 

k 2 ^ k 3 



on a Calabi-Yau manifold. In |Figure lj we fit = ff + ff for k > 3 and find good 



agreement with the data points. 

An important question is how many points are necessary to approximate the 
Calabi-Yau threefold in the numerical integration for any given k. The problem is that 
we really are trying to compute the Nk x iV^-matrix h a ^, whose dimension increases 
quickly with k, see ITable 1[ Hence, to have more equations than indeterminates, we 
expect to need 

N P > Nl (54) 
points to evaluate the integrand of the T-operator on. To numerically test this, we 
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k = k = l k = 2 k = 3 k = A k = 5 k = 6 k = 7 k = 8 

5 (k) 



Figure 2: The error measure for the balanced metric on the Fermat quintic 
as a function of k, computed by numerical integration with different 
numbers of points N p . In each case, we iterated the T-operator 10 
times and evaluated on 5,000 different test points. Note that we 
use a logarithmic scale for the <jfc axis. 
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Figure 3: The error measure o~k for the balanced metric on the Fermat quintic 
as a function of = number of entries in h a/3 € MatjVfcxJVj,- In 
other words, evaluating the T- operator requires N? scalar integrals. 
In each case, we iterated the T-operator 10 times and finally evalu- 
ated Ok using 5,000 different test points. We use a logarithmic scale 
for both axes. 
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compute cTfc using different numbers of points N p . The result is displayed in |Figure2 
where we used the more convenient logarithmic scale for a^. Clearly, the error measure 
a*, starts out decreasing with k. However, at some A^-dependent point it reaches a 



minimum and then starts to increase. In Figure 3 we plot the same a k as a function 



of Ni. This confirms our guess that we need N p > Ni points in order to accurately 



perform the numerical integration. One notes that the data points in Figure 2 seem to 
approach a straight line as we increase N p . This would suggest an exponential fall-off 

a k « 0.523e" a324fc . (55) 

It is possible, therefore, that the theoretical error estimate eq. (!53l) could be improved 
upon. 

So far, we have applied our procedure to the Fermat quintic Qp for simplicity. 
However, our formalism applies equally well to any quintic Q in the 101-dimensional 
complex structure moduli space with the proviso that, for a non-Fermat quintic, one 
must use the LC\Q method of choosing points. An important property of the programs 
that implement our procedure is that they make no assumptions about the form of the 
quintic polynomial eq. i^j. We proceed as follows. First, fix a quintic by randomly (in 
the usual flat distribution) choosing each coefficient C( n0inijn2in3in4 ) on the unit disk, 
see eq. (j2J). Then approximate the Calabi-Yau metric via Donaldson's algorithm and 



compute the error measure a k - In Figure 4, we present the results for a k for five 



randomly chosen quintics, and compare them to the Fermat quintic. We observe that 
the convergence to the Calabi-Yau metric does not strongly depend on the complex 
structure parameters. 



3 Group Actions and Invariants 

3.1 Quotients and Covering Spaces 

Thus far, we have restricted our formalism to quintic Calabi-Yau threefolds Q C P 4 . 
These are, by construction, simply connected. However, for applications in heterotic 
string theory we are particularly interested in non-simply connected Calabi-Yau man- 
ifolds where one can reduce the number of quark/lepton generations and turn on 
discrete Wilson lines [3H [32J, [331 EH [351 ESI [37]. Therefore, it is of obvious interest 
to compute the metrics in such cases. However, these manifolds are more compli- 
cated than hypersurfaces in projective spaces. In fact, any complete intersection in a 
smooth toric variety will be simply connected^- Therefore, we are usually forced to 
study non-simply connected Calabi-Yau threefolds Y, 

m(X) = n ^ l, (56) 

Note, however, that there are 16 cases of smooth, non-simply connected hypersurfaces in singular 
toric varieties [55] , 
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Figure 4 '• The error measure o~k as a function of k for five random quintics, 
as well as for the Fermat quintic. The random quintics are the sum 
over the 126 quintic monomials in 5 homogeneous variables with 
coefficients random on the unit disk. We use a logarithmic scale for 
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via their universal covering space Y and the free group action II : Y — > Y. 

In order to carry through Donaldson's algorithm on Y, we now need to generalize 
the notion of "homogeneous polynomials" to arbitrary varieties. As mentioned previ- 
ously, the homogeneous coordinates on the quintic QeF 4 can be interpreted as the 
basis of sections of the line bundle 0q(1), 

span{z , 21,22,23, 24} = H°(Q, 0q(1)). (57) 

The special property of 0q(1) is that it is "very ample", that is, its sections define an 
embedding 

$G 5 (i) : Q -> P 4 , x (-> [zo(x) : Zl (x) : z 2 (x) : z 3 (z) : z 4 (z)] . (58) 

Hence, we need to pick a "very ample" line bundle on Y in order to compute the 
metric there. Furthermore, to discuss Y, we will also need to "mod out" by the 
group action. It follows that the group II must act properly on the line bundle. In 
mathematical terms this is called an "equivariant line bundle" , and there is a one-to- 
one correspondence 



Il-equivariant 


n 


Line bundles 


line bundles on Y 


n* 


on Y. 



(59) 



Let us denote such a line bundle on Y by XL. We are specifically interested in the 
sections of this line bundle, since they generalize the homogeneous coordinates. The 
important observation here is that the sections of a LT-equivariant line bundle on Y 
themselves form a representation of II. Furthermore, the Il-invariant sections corre- 
spond to the sections on the quotient. That is, 

H°(Y,L) U = H°(Y,L/U). (60) 

Hence, in order to compute the metric on the quotient Y = Y/Tl, we can work on the 
covering space Y if we simply replace all sections by the n-invariant sections. 

In this paper, we will always consider the case where Y is a hypersurface or a 
complete intersection in (products of) projective spaces. Then 

• The sections on the ambient projective space are homogeneous polynomials. 

• The sections on Y are the quotient of these polynomials by the defining equa- 
tions. 

• The invariant sections on Y are the invariant homogeneous polynomials modulo 
the invariant polynomials generated by the defining equations. 

The mathematical framework for counting and finding these invariants is provided by 
Invariant Theory [22], which we review in the remainder of this section. 
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3.2 Poincare and Molien 

Let C[x] be a polynomial ring in n commuting variables 

x= (xi, . . . ,x n ). (61) 

As a vector space over the ground field C, it is generated by all monomials 

C[x] = CI © Cxi © ■ • • © Cx n © Cxi © ■ ■ • • (62) 

Clearly, C[x] is an infinite dimensional vector space. However, at each degree k we have 
a finite dimensional vector space of homogeneous degree-k polynomials. A concrete 
basis for the degree-/c polynomials would be all distinct monomials of that degree. 

By definition, the Poincare series is the generating function for the dimensions of 
the vector subspaces of fixed degree, that is, 

oo 

P{C[x],t) = (dim c C[f] fc ) t k (63) 

k=0 



where C[x\k is the vector subspace of C[x] 
mial ring in n commuting variables x±, . . 
bosons, and one obtains 

n 



of degree k. The monomials of the polyno- 
, . , x n can be counted just like n species of 

k=0 ^ ' 



We have already mentioned that the homogeneous degree-fc polynomials in n variables 
are just the sections of P n-i(/c). Hence, the number of degree-/;; polynomials is the 
same as the dimension of the space of sections of the line bundle P n-i(A;), 

dim c C[f] fc = /i°(P n ~ 1 ,0 P n-i(A;)). (65) 

Acknowledging this geometric interpretation, we also write 

oo 

P((V-i,t) = ^/i°(P n - 1 ,0 P n-i(A;))t fc = P(C[x\,t). (66) 

k=0 

Furthermore, note that 

P(M © M',t) = P(M,t) + P(M',t) (67) 

for any rings M and M' . 

A n-dimensional representation of a finite group G generates a group action on the 
polynomials eq. (IB"!?]) . One is often interested in the invariant polynomials under this 
group action, which again form a ring C [x\ G . Clearly, the invariant ring is a subring 
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of C [x]. Since the group action preserves the degree of a polynomial, one can again 
define the Poincare series of the invariant ring, 

oo 

P{C[xf,t) =^(dim c C[f]^)t fc . (68) 

k=0 

The coefficients in eq. (|68|) can be obtained using 

Theorem 3 (Molien). Let G C GL(n, C) be a finite matrix group acting linearly on 
the n variables x = (x%, . . . ,x n ). Then the Poincare series of the ring of invariant 
polynomials, that is, the generating function for the number of invariant polynomials 
of each degree, is given by 

1 1 geG v a > 

Equation (169!) is called the Molien formula. 



3.3 Hironaka Decomposition 

Although eq. fl69|) contains important information about C[xi, . . . , x n ] G , the most de- 
tailed description is provided by the Hironaka decomposition, which we discus next. 
To construct this, one first needs to find n homogeneous polynomials 9\, . . . , 9 n , in- 
variant under the group action, such that the quotient 

C[x 1 ,...,x n ]/(9 1 ,...,e n ) (70) 

is zero-dimensional. The above condition is equivalent [3H] to demanding that the 
system 9i — 0, i — 1, . . . n has only the trivial solution. This guaranties that the 9i 
are algebraically independent. Then 

Theorem 4 (Hironaka decomposition). With respect to 9\, . . . , 9 n chosen as above, 
the ring of G-invariant polynomials can be decomposed as 

C[xf = 77iC[0i, . . . A] © V2C[9i, • • • , 9 n ] © • • • © r, s C[9x, . . . , 9 n \. (71) 

Clearly, the r]i are themselves G-invariant polynomials in C[x]. Thus any G- 
invariant polynomial is a unique linear combination of rj^s, where the coefficients 
are polynomials in #j. The polynomials 0, are called the "primary" invariants and 
r)j the "secondary" invariants. Note that, while the number of primary invariants is 
fixed by the number of variables x\, . . . , x n , the number s of secondary polynomials 
depends on our choice of primary invariants. Using eq. fl64l) with each Xi replaced by 
9i, we find that the Poincare series for C[9\, . . . , 9 n ] is given by 

P(C[0i, . . . , 9 n ],t) = —- td 1 tdcg(en)) • (72) 
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Moreover, multiplication by rji shifts all degrees by deg(?7j). Therefore, applying 
eq. (1671) we obtain the Poincare series for the Hironaka decomposition, 



t D, t D a 



'1 _ t dl) . . . (1 - ^n) (1 _ fdl) . . . (I _ t dn 

t D1 + • • • + 1 



Dl , . , +Ds (73) 



'1 - t d A ■■■(!- t d - 



where Dj = deg(rjj) and d{ = deg(^). Each term in the numerator of eq. (1731) 
corresponds to a secondary invariant. 



4 Four-Generation Quotient of Quintics 
4.1 Four Generation Models 

Were one to compactify the heterotic string on a generic quintic Q using the standard 
embedding, then the four-dimensional effective theory would contain \x{Q) — 100 
net generations. A well known way to reduce this number [I] is to compactify on 
quintics that admit a fixed point free Z 5 x Z 5 action. In that case, the quotient 
manifold Q = Q j (Z 5 x Z 5 ) has only \x(Q) — [zixld = ^ g enera tions. In this section, 
these special quintics and their Z 5 x Z 5 quotient will be described. We then compute 
the Calabi-Yau metrics directly on these quotients Q using a generalization of our 
previous formalism. 

Recall from ISection 2l that a generic quintic Q C P 4 is defined as the zero locus of 
a degree-5 polynomial of the form eq. ([2]). In general, it is the sum of 126 degree-5 
monomials, leading to 126 coefficients C( no nijn2in3jn4 ) G C. However, not all of these 
quintic threefolds admit a fixed point free Z 5 x Z 5 action. To be explicit, we will 
consider the following two actions on the five homogeneous variables defining P 4 , 
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(74) 



Clearly g\ — 1 = g%, but they do not quite commute: 

9i92 = e 2 f L g 2 gi Q^gT 1 g^ 1 = ■ (75) 
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However, even though g\ and g% do not form a matrix representation of Z5 x Z5, they 
do generate a Z 5 x Z 5 action on IP because on the level of homogeneous coordinates 
we have to identify 



2iri 2iri 2ivi 2ni 2iri 



(77) 



[z Q : z 1 : z 2 : z 3 : z 4 ] = [e z z : e * z 1 : e z 2 : e z 3 : e 3 z 4 ] 
= 91929Y 1 92 1 ([z Q : zi : z 2 : z 3 : zM. 

If the quintic polynomial Q{z) is Z 5 x Z 5 -invariant, then the corresponding hypersur- 
face will inherit this group action. One can easily verify that the dimension of the 
space of invariant homogeneous degree-5 polynomials is 6, as we will prove in eq. fl9Tj) 
below. Taking into account that one can always multiply the defining equation by a 
constant, there are 5 independent parameters 0i, . . . , 05 G C. Thus the Z 5 x Z 5 sym- 
metric quintics form a five parameter family which, at a generic point in the moduli 
space, can be written as 

Q{z) = (zl + z\ + z\ + z\ + zl) 
+ 01 (z zi z 2 z 3 z 4 ) 

+ 02(^0^1^4 + Zqz\z 2 + z z 3 zl + Z\z\z 3 + z 2 z\z^ 

,i/2 2,2 2,2 2,2 2,2 2\ 
t ( P3\ z o z l z 2 ' z i z 2 z 3 ~r Z 2 Z 3 Z 4 : + z 3 z i z ' z 4 z z l) 

,i(2 2 ,22 ,22 ,22 ,22\ 
T" (P4\ z o z l z 3 1 Z \ Z 2 Z A 1 Z 2 Z 3 Z ' Z 3 Z 4 Z 1 + z i Z Q Z 2) 

+ 05(^0^2^3 + z i z 3 z 4: + Z 2 Z 4 Z + Z 3 Z Z 1 + z \ z l z 2)- 

The explicit form of these invariant polynomials is derived in lSubsection 4.3l and given 
in eq. (I9T1) . Note that, even though the Z 5 x Z 5 action on P necessarily has fixed 
points, one can check that a generic (that is, for generic 0i, . . . , 0s) quintic threefold 
Q is fixed-point free. 

Now choose any quintic defined by eq. ( 1771) . Since the Z 5 x Z 5 action on it is fixed 
point free, the quotient 

Q = Q/(Z 5 xZ 5 ) (78) 
is again a smooth Calabi-Yau threefold. Its Hodge diamond is given by [40] 

1 



/- .010 
h?>*(Q) =/Mq/(Z 5 xZ 5 )) = 1 5 5 1, (79) 

1 


1 

where we again see that there is a h 2,1 (Q) = 5-dimensional complex structure moduli 
space parametrized by the coefficients 0i, . . . , 05. 
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4.2 Sections on the Quotient 

We now extend Donaldson's algorithm to compute the Calabi-Yau metric directly on 
the quotient Q = Q j (Z 5 x Z 5 ). To do this, we will need to count and then explicitly 
construct the Z5 x Z5 invariant sections, that is, the Z5 x Z5 invariant polynomials, on 
the covering space Q G P 4 , as discussed in ISubsection 3.11 These then descend to the 
quotient Q and can be used to parametrize the Kahler potential and the approximating 
balanced metrics. 

One technical problem, however, is that the two group generators g\ and g 2 in 
eq. (|74|) do not commute; they only commute up to a phase. Therefore, the homoge- 
neous coordinates 

span {so.jzi,^, a*, *4 =H°(Q,Qq(1)) (80) 

do not carry a Z5 x Z5 representation. The solution to this problem is to enlarge the 
group. Each generator has order 5 and, even though they do not quite generate Z 5 x Z 5 , 
they commute up to a phase. Hence, g\ and g 2 generate the "central extension" 

1 — ► Z 5 — ► G — > Z 5 x Z 5 — ► 1 (81) 

with |G| = 125 elements. This group G is also called a Heisenberg group since it is 
formally the same as [x,p] = 1, only in this case over Z 5 . It follows that H°(Q, CWl)) 

does carry a representation of G and, hence, so does H°(Q, Og(k)) for any integer k. 

Note that, when acting on degree-k polynomials Pk(z), the commutant eq. (!75|) 
becomes 

9i929T 1 92 1 {pk{z)) = e 2 ^ p k (z). (82) 

Therefore, if and only if k is divisible by 5 then the G representation reduces to a true 
Z 5 x Z 5 representation on H°(Q, Oq{k)). That is, k must be of the form 

k = 5£, £eZ. (83) 

The formal reason for this is that only the line bundles 0^(5^) are Z 5 x Z 5 equivariant. 
The invariant subspaces of these Z 5 x Z 5 representations define the invariant sections. 
Hence, we only consider homogeneous polynomials of degrees divisible by 5 which are 
invariant under the action of Z5 x Z5 in the following. 

4.3 Invariant Polynomials 

As a first step, determine the Z5 x Z5 invariant sections on the ambient space P . 
That is, we must find the invariant ring 

C[zq, 3l, 3 2 , 3 3 , 3 4 ] G (84) 
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over P 4 , where G is the Heisenberg group defined in the previous subsection. One can 
read off the number of invariants Njf at each degree k from the Molien series 



P{C[z , Zl , z 2 , z 3 , z,] G , t)=J2^t k = T^Yl 

k ' ' geG 



det (l - tg) 



= 1 + 6t 5 + 41t 10 + 156t 15 + 426t 20 + 951t 25 + 1856t 30 + 3291t 35 + 5431t 40 + 

+ 8476t 45 + 12651t 50 + 18206t 55 + 25416t 60 + 34581t 65 + • • • (85) 

We see that the only invariants are of degree k = 5£, as discussed in the previous 
subsection. To go further than just counting the invariants, one uses the Hironaka 
decomposition which was introduced in ISubsection 3.3l For that, we need to choose 
5 primary invariants, the same number as homogeneous coordinates. Unfortunately, 
any 5 out of the 6 quintic invariant polynomials are never algebraically independent. 
Hence, picking five degree-5 invariants never satisfies the requirements for them to be 
primary invariants. It turns out that the primary invariants of minimal degree consist 
of three degree-5 and two degree-10 invariants, which we will list in eq. (1881) below. 
First, however, let us rewrite the Molien series as in eq. (173D . 



, r \ 1 + 3t 5 + 24t 10 + 44t 15 + 24t 20 + 3t 25 + t 30 , . 

P(C[z ,z 1 ,z 2 ,z 3 ,z i \ ,t) = -3- -3 . (86) 

; (1 - t 5 ) 3 (l - t 10 ) 

We see that this choice of primary invariants requires 

1 5 ^ 3 10 2 
— TTdegfli = —— = 100 = 1 + 3 + 24 + 44 + 24 + 3 + 1 (87) 

secondary invariants in degrees up to 30. We again note that this decomposition is not 
unique, as one can always find different primary and secondary invariants. However, 
our choice of primary invariants is minimal, that is, leads to the least possible number 
(= 100) of secondary invariants. 

Knowing the number of secondary invariants is not enough, however, and we need 



the actual G- invariant polynomials. As will be explicitly checked in Appendix A the 
five G-invariant polynomials 

61= zl + z\ + zl + zl + zl =4 + ( C Y C ) 

62 = ZqZxZ^Z-^Z^ 

6 3 = zlz x z A + z zfz 2 + z z 3 zl + Z\z\z 3 + z 2 z\z± = zlz x z A + (eye) (88) 

6 4 = Z™ + zf + zf + zf + zf = zl° + (eye) 

O5 = zlz 2 z 3 + ZqZxzI + 20^4 + z\z 3 z A + z x z 2 z\ = z\z 2 z 3 + (eye) 

satisfy the necessary criterion to be our primary invariants, where (eye) denotes the 
sum over the five different cyclic permutations z — > z\ —> ■ ■ ■ — > Z4 — > z . Next, 
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we need a basis for the corresponding secondary invariants, which must be of degrees 
0, 5, 10, 15, 20, 25, and 30 according to eq. (1861) . In practice, these 100 secondary- 
invariants can easily be found using Singular |^Tl They are 



T] 2 = z\z x z\ + (eye), r] 3 = zlz\z % + (eye), 7/4 = z\z 2 z z + (eye), 



775 = 44 + (cyc), 
vs = 444 + ( c y c ); 

Til = zlz\z\ + (eye), 
rj 14 = z Q z\z 2 + (eye), 
7717 = z\z\ z 4 + (eye), 

V2Q = 4444 + ( c yc), 
V23 = ZQzfzjz 3 + (eye), 
V26 = 444 z * + (eye), 

??29 = zf + (eye), 
V32 = ZQ°zf + (eye), 
7/35 = zlz\z\ + (eye), 
(eye), 

(eye), 



^38 — 444 



44 z 3 + (eye), 



8 4 

r)4i = z z 2 z_ 

r/ 47 = 444 + (eye) 
r/50 = zllzxzl + (eye), 
7753 = z\lz x z\ + (eye), 

Vsb = 4444 + (eye) 
r/59 = zlz\z\z\ + (eye) 
%2 = 444 z * + (eye) 
??65 = 4444 + (eye) 
Tes = 444 z ± + (eye) 
7771 = zlzizlz A + (eye) 



'/« 
'/!) 
7/12 
Tl5 
?7l8 
?721 
7?24 
7?27 

??30 
??33 
7/36 
7/39 
7?42 
7/45 
??48 
7?51 
7?54 
^757 
??60 
??63 
7/66 
7/69 
7/72 



z o44 + 

z A z 4 z 2 + 

^0^1-^3 ^ 
^0^1 Z 3 + 

^O 2 ^ 2 ^ + 

7 2 1 

^0 Z 1 2 3 + 

^0 Z l Z 3 Z/ i 

4 4 
Zq ZiZ 2 Z 3 

4z\z\z 3 



(eye), 
(eye), 
(eye), 
(eye), 
(eye), 
+ (eye), 7722 
+ (eye), 7725 
+ (eye), 7728 



4° 4 + (eye), 
444 + (eye), 

7 4 4 
Z Q Z x £ 4 



r 8 r r 6 



y9 --/4 2 

'0^1^3 
,11 ^3, 



(eye), 
444 + (eye), 

*o*i*i + (cyc), 

(eye), 
(eye), 
Zq1z(z 2 + (eye), 
z\2z\z z + (eye), 

4444 + (eye) 

4444 + (eye) 
444 z 3 + (eye) 
zlz\z 2 z\ + (eye) 

Zo^i ^3 + (eye) 

zo«i^4 + ( c y c ) 



7/7 
7/10 
7/13 
7?16 
7/19 



7/31 
7/34 
7/37 
7/40 
7/43 
7/46 
7/49 
7/52 
7/55 
7/58 
7/61 
7/64 
7/67 
7/70 



z o44 + 
444 + 

Z Z 2 Z 4: 
7 2 1 

ZqZ\Z4 + 
^0 ^1 ^2 ^4 
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„ _ JO, / \ _ r 30 1 -30 1 -30 , 30 , „30 

^100 — z o + icycj — z + z x + z 2 + z 3 + z 4 . 

The Hironaka decomposition of the ring of G-invariant homogeneous polynomials is 
then 
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C[zq, Zi,Z 2 , 2 3 , 2 4 ] G = C[0i,0 2j 03i 04,0s]- 



(90) 



i=l 



As a simple application, we can read off a basis for the invariant degree-5 polynomials, 

(91) 



C[z , 21, 22,23,24] 5 = span |?7i0i, 7716*2, r/i0 3 , 772, ?73, ?74 
= span 1 1; 2 , 3 , 772, ?73, tytj- 



Note that this is the basis of invariant quintic polynomials used in eq. (ITT)) to define 
Q(z). 

4.4 Invariant Sections on the Quintic 

The next step is to restrict the G- invariant sections on P 4 to the hyper surf ace Q. 
In ISection 2| we showed how to accomplish this for all sections on generic quintics 
Q G P 4 . Since the sections on the ambient space are nothing but homogeneous 
polynomials, the restricted sections were the quotient of the homogeneous polynomials 
by the hypersurface equation Q = 0, 



#°(p 4 ,(V(fc)) 



C[z , 2i, 22, 23, 24J 



restrict 



H°(Q,Odk)) 



(92) 



C[2 , 21,22,23,24] / (Q) 



Now consider the quintics defined by eq. (1771) . which allow a Z 5 x Z 5 action. Here, 
one only wants to know the G-invariant sections on Q, since these correspond to the 
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sections on the Z5 x Z5 quotient Q = Q j (Z5 x Z5). Moreover, since the G-invariant 
polynomials are of degree 5£, we only consider this case. Hence, the G-invariant 
sections are 



#°(P 4 ,(V(5£)) 



G 



restrict 



H°(Q, 0q(5£)) 



G 



(93) 



C[z 



[C[Zo,Z 1 ,Z 2 ,Z 3 ,Z 4 \/ (Q)) 



Finally, we identify the invariant sections on Q with sections on the quotient manifold 
Q, as discussed in ISubsection 3.11 Therefore, the sections on Q are 



H°[Q,Oq(5£)/(Z 5 xZ 5 )\ =H\Q,Qq{W)Y = [C^zuz^z^zt] {Q)\ . (94) 



By unravelling the definitions and using eq. (190]) . the invariant subspace of the quotient 
ring is given by 



C[z ,z 1 ,z 2 ,z 3 ,z 4 ] / (Q(z)) 



C[z , Zi, z 2 , z 3 , z A } G (Q(z)) 



G 



100 



100 



07/iC[0i,0 2 ,03,04,0 5 ] / 0Q^C[0!, 2 , 03,04,05] ■ (95) 



i=l 



i=l 



Using eq. (ITT)) , the hypersurface equation is 

Q{z) = 

zl + z\ + zl + z\ + z\= - 4>i [zqZ\Z 2 z 3 z±) - 



& (96) 



01 = - 0102 - 0203 - 03% - 04% ~ 05%, 



and, hence, we can simply eliminate d\. Therefore, forming the quotient is particularly 
easy, and we obtain 



g 100 

C[z , *1, *2, Z 3 , Z A ] I (Q(Z))) = Tji C[0 2 , 3 , 4 , 5 ]- 



(97) 



i=i 



We list the number iv"^ of G-invariant degree-5£ polynomials on P 4 as well as the 
number of invariant polynomials after restricting to Q, N^ e , in ITable 21 Since we 
know the homogeneous degrees of the primary and secondary invariants, and r\ 
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u 


5 


10 


15 


20 


25 


30 


35 


40 


Ng 


6 

5 


41 

35 


156 
115 


426 
270 


951 
525 


1856 
905 


3291 
1435 


5431 
2140 



Table 2: The number of G -invariant degree hi -homogeneous polynomials Ng, 
eq. (|85p , and the number of remaining invariant polynomials Ng after 
imposing the hypersurface equation Q(z) = 0, see eq. (fT7|) . 



respectively, it is a simple combinatorial problem to list all monomials in eq. ([9 
of fixed degree 5£. They then form a basis for the sections on Q, 



N G -1 



H°(q, Oq(5£)/(Z 5 x Z B )) = Bpan{s a }2 

f 100 \ 100 

^ViC[9 2 ,9 3 ,9 4 ,9 5 \) =Q)viC[92,9 3 ,9 4 ,9 5 ]M_ degTli . (98) 



1=1 / « i=l 



4.5 Results 

We have now computed an explicit basis of invariant sections of Qq(5£), which can be 

identified with a basis of sections on the quotient manifold Q = Q/(Zs x Z5). This 
is all we need to extend Donaldson's algorithm to Q. Literally the only difference in 
the computer program used in ISubsection 2.41 is that now 

• the degree of the polynomials must be k = 5£, £ G Z>, and 

• the sections are given in eq. (1981) . 

Hence, one can compute the balanced metrics on Q. As £ — > 00, these will approach 
the unique Calabi-Yau metric. We write a$z(Q) f° r the error measure computed 
directly for the balanced metrics on the non-simply connected threefold Q. Note 
that there is still a 5-dimensional complex structure moduli space of such threefolds. 



However, as we have seen in Figure 4 , the details of the complex structure essentially 
play no role in how fast the balanced metrics converge to the Calabi-Yau metric. 
Therefore, as an example, in |Figure~5] we plot for the quotient Qf = Qf/ (Z5 x Z 5 ) 



of the Fermat quintic. Note that the error measure tends to zero as £ — > 00, as it 
should. 

Comparison With the Covering Space 

We have now extended Donaldson's algorithm so as to compute the successive approx- 
imations to the Calabi-Yau metric directly on the quotient manifold Q. Clearly, these 
metrics can be pulled back to Z 5 x Z 5 symmetric metrics on the covering space Q, 
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Figure 5: The error measure cr^Qp) on the non-simply connected threefold 
Qf = Qf/ (^5 x ^5)- For each £ € TL > we iterated the T-operator 10 
times, numerically integrating using N p = 1,000,000 points. Then 
we evaluated o~m{Qf) using 20,000 different test points. Note that 
all three plots show the same data, but with different combinations 
of linear and logarithmic axes. 
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k 

5 10 15 20 25 




5 10 15 20 25 150,000 300,000 

U N 2 = # of entries in h a ^ e Mat^*^ 

Figure 6: The metric pulled back from Qf = Qf/{^5 x ^5) compared with 
the metric computation on Qf- The error measures are ct^Qf) 
and <ta:(Qf); respectively. On the left, we plot them by the degree 
of the homogeneous polynomials. On the right, we plot them as a 
function of N 2 , the number of sections squared. On Qf, the number 
of sections is N^; on Qf the number of sections is N^. The a -axis 
is logarithmic. 



32 



thus approximating the Calabi-Yau metric on Q. Let us denote by iou the balanced 
Kahler form on Q computed at degree-5£, and by q*u 5 e its pull-back to Q. We define 
&u(Q) to be the error measure evaluated using the pull-back metric, that is, 



VoIcy (Q) Jq 



1 - 



<Z*u&/Vol K (Q) 



QAQ/ VoIcy (Q) 



dVol 



CY 



VoIcy(Q) 



Q 



1 - 



Vol K (Q) 



o a n / voIcy(Q) 



dVol 



CY 



(99) 



Now recall that in ISection 21 it was shown how to determine the Calabi-Yau metric 
on any quintic threefold. This, of course, includes the Z 5 x Z 5 quintics Q defined 
by eq. ( 1771) . However, since most quintics do not admit a finite group action, the 
procedure specified in ISection 21 finds the Calabi-Yau metric using generic homoge- 
neous polynomials. That is, it finds an explicit polynomial basis for H°(Q,Og(k)), 
computes the balanced metric and determines the Calabi-Yau metric as the k — > oo 
limit. When applied to our Z 5 x Z 5 quintics, this second method will also compute 
the unique Z 5 x Z 5 symmetric Calabi-Yau metric. However, it does so as the limit 
of balanced metrics constructed from sections of Og(fc) which do not share this sym- 
metry, rather than from invariant sections of 0q(5£) as above. That is, this second 
method does not exploit the Z5 x Z5 symmetry. The associated error measure <Tfc 
is evaluated using eq. ( 1491) for Z5 x Zs-symmetric quintics. It is of some interest to 
compare these two methods for calculating the Calabi-Yau metric on Q. Specifically, 



in the left plot of |Figure 6] we compare the error measure to on the Fermat 
quintic Qp. Interestingly, for fixed degrees k = 5£ the pull-back metric is a worse ap- 
proximation to the Calabi-Yau metric on Q than the metric computed on Q without 
taking the symmetry into account. The reason is that, in addition to the Z 5 x Z 5 
invariant polynomials on Q, there are many more that transform with some character 
of Z 5 x Z 5 . These polynomials provide extra degrees of freedom at fixed degree 5£, 
which allow the balanced metric to be a better fit to the Calabi-Yau metric. 

However, a more just comparison is by the amount of the numerical effort, that is, 
the number (iV^) 2 and (A^) 2 , respectively, of entries in the h a/3 matrix. We plot 
and <7fc as a function of N 2 



cases 5£ = k = 5 and 5£ 



m 

k 



Figure 6 We see that, except for the two lowest-degree 



10, the pull-back metric computation (that is, using 



invariant sections) is more efficient. 
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5 Schoen Threefolds 



5.1 As Complete Intersections 

By definition, Schoen type Calabi-Yau threefolds are the fiber product of two dPg 
surfaces, B 1 and B 2 , fibered over P 1 . Recall that a dPg surface is defined as a blow-up 
of P 2 at 9 points. In principle, these points can be "infinitesimally close", that is, 
one of the blow-up points lies within a previous blow-up, but we will only consider 
the generic case where all 9 points are distinct. Moreover, we are going to restrict 
ourselves to the case where "no Kodaira fibers collide". In that case, the Hodge 
diamond of the Schoen threefold X is [U [36j S3] 

1 


19 

h p ' q {X) =1 19 19 1 . (100) 
19 

1 

These generic Schoen Calabi-Yau threefolds can be written as a complete intersec- 
tion as follow [HI [131 [33 IB]. First, consider the ambient variety P 2 x P 1 x P 2 with 
coordinates 

([x : Xl : x 2 ], [t : t x ), [y : Vl : y 2 ]j E P 2 x P 1 x P 2 . (101) 

The Calabi-Yau threefold X is then cut out as the zero-set of two equations of multi- 
degrees (3, 1, 0) and (0, 1,3), respectively. The two equations are of the form 

P(x,t,y) = t Pi (x Q ,xi, x 2 ) + hP 2 (x ,x 1 ,x 2 ) = 0, (102a) 
R(x, t, y) = ti-Ri (y , yi, y2) + hFL 2 (y , y h y 2 ) = (102b) 

where Pi, P 2 , Ri, and R 2 are cubic polynomials. The ambient space P 2 x P 1 x P 2 is a 
toric variety and X is a toric complete intersection Calabi-Yau threefold [ISl HHJ HTj. 

5.2 Line Bundles and Sections 

The first Chern classes of line bundles on X form a 

h 1 ' 1 (X) = 19 (103) 

dimensional lattice. Note, however, that most of them do not come from the ambient 
space which has 

/ i 1 ' 1 (P 2 xP 1 xP 2 ) =3. (104) 
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In other words, most of the divisors D and their associated line bundles -C(-D) are not 
toric; that is, they cannot be described by toric methods. We could embed X in a 
much more complicated toric variety [16] where all divisors are toric. However, for 
now^) we will simply ignore the non-toric divisors and restrict ourselves to line bundles 
on X that are induced from P 2 x P 1 x P 2 . 

The line bundles on P 2 x P 1 x P 2 are classified by their first Chern class 

ci^O P 2 xP i xp2 (a 1 ,6,a 2 )^ = (ai,6,a 2 ) e Z 3 = H 2 (F 2 x P 1 x P 2 ,Z). (105) 

Just as in the P 4 case previously, their sections are homogeneous polynomials of the 
homogeneous coordinates. Now, however, there are three independent degrees, one for 
each factor. That is, the sections of P 2 x p i x P 2 (ai, b, a 2 ) are homogeneous polynomials 
of 

• degree a\ in xo, x±, x 2 , 

• degree b in t , t\, 

• degree a 2 in y , y ly y 2 . 

The number of such polynomials (that is, the dimension of the linear space of poly- 
nomials) is counted by the Poincare series 

p(0 P 2 x¥ i x¥ 2,(x,t,y)j = ^°(P 2 xP 1 xP 2 ,0(a 1 ,6,a 2 ))x ai tV 2 

ai ,b,a2 (106) 
111 v ; 



(i-x)3(i-t)2(i-y) 8 ' 



We now want to restrict the sections to the complete intersection X C P 2 x P 1 x P 2 ; 
that is, find the image 

H° ( P 2 x P 1 x P 2 , 0(a 1; b, a 2 )) rcstrict > H° (x, 0$ (a u b, a 2 )) (107) 

foi@ ai,6, a 2 > 0. As discussed previously, this amounts to finding a basis for the 
quotient space 

H°(X,Qx(ai,b,a2)) = (c[x Q , x u x 2 ,t , h, y , y h y 2 ] /(?, R)) (108) 

\ / / (ai ,b,a,2) 



8 This will be partially justified in IScction"6l where we investigate a certain Z3 x Z3-quotient of 
X. There, only the toric line bundles will be relevant. 

9 Note that ci (Oj^(ai, b, 02)) 6 if 2 (X,Z) is in the interior of the Kahler cone if and only if 
01, b, a-i > 0, see H51 . 
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of degree (01,6,02). Note that this quotient by more than one polynomial is much 
more difficult than the case where one quotients out a single polynomial, as we did for 
quintics in ISection 21 In general, this requires the technology of Grobner bases [39] . 
Suffices to say that we are in a very advantageous position here. 

By a suitable coordinate change, we can assume that the t 2/o term in R is absent. 
That is, 

P = t X 3 n H 

~ 3 2 (109) 

R=0- t y* + hVoVi + ■■■ ■ 

Then, for otherwise generic polynomials P and R and lexicographic monomial order 

xo -< yo -< to -< xi -< yi -< h -< x 2 -< y 2 , (HO) 

the two polynomials generating 

(P,R) C C[x ,xi,X2,t ,ti,y ,y 1 ,y 2 \ (111) 

already form a Grobner basis. This means that the quotient in eq. ( 1108j) can be 
implemented simply by eliminating the leading monomials to^o an d toVoVi i n the 
polynomial ring C[x 0l x 1 ,X2,to,t 1 ,y ,y 1 ,y 2 \. 



5.3 The Calabi-Yau Volume Form 

As in the case of a hypersurface, one can express the (3, 0)-form of the complete 
intersection as a Griffiths residue. By definition, the zero loci P = and R = 
intersect transversally, so one can encircle each in an independent transverse direction. 
The double residue integral 

^ff(fxdt(fy 

^ff^nr < 112 > 

is again independent of the chosen inhomogeneous coordinate chart. Hence, it defines 
a holomorphic (3, 0)-form which must be the holomorphic volume form. 



5.4 Generating Points 

Since the defining Equations (1102aj) . (1102bl) are at most cubic in the x and y coordi- 
nates, there is a particularly nice way to pick points. This is a generalization of the 
L H Q method presented in ISubsection 2.31 to generate points in generic quintics. In 
the present case, select a specific P 1 x P 1 in the ambient space, namely, 

P^pt.jxP 1 c P 2 xP : xP 2 . (113) 
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This can easily be done with an 577(3) x SU(2) x SU (3) -invariant probability density 
of such configurations. The intersection 



(V x{pt.} x P 1 ) HI = {9 points} (114) 

consists of nine points. To compute the coordinates of the nine points, one needs to 
solve two cubic equations, which can be done analyticallj@ 

We still need the distribution of these "random" points. First, note that there are 
three obvious (1, l)-forms. These are the pull-backs 

7r*(^ p2 ), 7r;(w P i), 7r;(cj p2 ) (115) 

of the standard (SU(m + 1) symmetric) Fubini-Study Kahler forms on P m , where 
Hi is the projection on the z-th factor of the ambient space. However, here the 
S77(3) x SU{2) x SU(3) symmetry of the ambient space is not enough to determine 
the distribution of points uniquely. 

In general, the question about the distribution of zeros was answered by Shifman 
and Zelditch [50]. Let us quickly review the result. Let £ be a line bundle on a 
complex manifold Y and pick a basis sq, . . . , sn-i of sections 

span{ So ,...,SAr-i} = #°(y,£). (116) 

Moreover, let £ be base-point free, that is, the sections do not have a common zero. 
In other words, 

$ £ : Y - P^" 1 , x i ► [s (t) : Sl (t) : SjV -i(*)] (117) 

is a well-defined map. The sections generate the TV-dimensional vector space H°(Y, £>) 
which contains the unit sphere SH°(Y, £). In other words, if we define sq, . . . , sn-i 
to be an orthonormal basis, then SH°(Y, L) is the common SU(N)-orbit of the basis 
sections. We take a random section s G SH°(Y, £) to be uniformly distributed with 
respect to the usual "round" measure, that is, SU (TV)-uniformly distributed. 

Finally, switch from each such section s to its zero locus Z s in Y, and consider the 
expected distribution of the random zero loci. Then 

Theorem 5 (Shifman, Zelditch). Under the above assumptions (in particular, that 
$£ is well-defined) the expected distribution of zero loci Z s is 

(Zc) = ^I^fs, (118) 



where ujps is the standard Fubini-Study Kahler form on P 



N-l 



10 Rccall that, to generate points on the quintic, we had to solve a quintic polynomial. This can 
only be done numerically. 
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(ai,b, a 2 ) 


(1,1,1) (2,2,2) 


(3,3,3) 


(4,4,4) 


(5,5,5) 


(6,6,6) 


N(ai ,6,02) 


18 108 


400 


1125 


2646 


5488 




18 108 


343 


801 


1566 


2728 



Table 3: The number of degree (01,6,02) -homogeneous polynomials N^ ai b ct2 ^ 
over P 2 x P 1 x P 2 and the number of remaining polynomials Nr ai ^ a2 ^ 

on X after imposing the two equalities P = = R defining the com- 
plete intersection. 

Note that, in our case, the embedding X C P 2 x P 1 x P 2 is generated by the three 
line bundles 

i/°(X,O^(l,0,0))=span{x ,x 1 ,x 2 } $ 0jf (1,0,0) : X P 2 , 

if°(X,O^(0,l,0))=span{t ,ti} $0,(0,1,0) : X - P 1 , (119) 

H°(X, 0^(0,0, l))=span{y , 2/1,2/2} $0,(0,0,1) : X ^ P 2 . 

Although none of the three $ maps is an embedding, they are all well-defined. This 
is sufficient for the theorem of Shifman and Zelditch. We point out that the $ maps 
are nothing but the restriction of the projections it to X C P 2 x P 1 x P 2 , 

$Ojf (1,0,0) = TTlIx: $0,(0,1,0) = ^Ijp $0,(0,0,1) = ^3 lx- (120) 

Hence, the expected distribution of a zero-loci of sections on X is 

(^0,(1,0,0)) ~ ^(^pOIx' (^Ox(o,i,o)) ~ ^(^pOIx' (^%(oai)) ~ ^(^Oljr 

(121) 

These are precisely the three (1, l)-forms we introduced previously in eq. f 1 1 1 5 j) . There- 
fore, if we independently pick the two P 1 factors and the point in eq. (II 13j) . then the 
distribution of simultaneous zero loci is 

dA ~ 7T* (cj p2 ) A vr^cjpi) A ^3(^2) (122) 

In other words, the points generated by the above algorithm are randomly distributed 
with respect to the auxiliary measure dA 



5.5 Results 

The new feature of the Schoen Calabi-Yau threefold, as opposed to the quintic, is that 
one now has different directions in the Kahler moduli space. On quintic threefolds 
there is only one Kahler modulus, which is just the overall volume. Now, however, 
there is a 19 = h l,1 (X) dimensional Kahler moduli space of which we parametrize 
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3 directions by the toric line bundles 0^(01,6,02). Note that, here as elsewhere in 
algebraic geometry, one has to work with integral Kahler classes that are the first 
Chern classes of some line bundle. This is not a real restriction, however, since any 
irrational slope direction in the Kahler moduli space can be approximated by a rational 
slope. A line with rational slope always intersects points in if 2 (X,Z). 



0.1 



0.01 





1 1 

AT _ 1 1\t2 _ 
1 JV (1,1,1) — io ; iV (l,l,l) ~~ 


1 1 

&(k,k,k) ' 1 ' 
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\ ^(2,2,2) = 


108, Nf 2 ^ 2) = 11,664 




^(3,3,3) 


= 343, iVjgg) = 117,649 




1 1 


1 1 



A; = k=l k = 2 k = 3 k = A k = 5 

Ox(k,k,k) 



Figure 7: The error measure OYfc^fc) for the metric on a Z3 x Z3 Schoen three- 
fold X . We iterated the T-operator 5 times, numerically integrating 
using N p = 1,000,000 points. Finally, we integrated <T(fc,fc,fc) using 
10,000 points. -/V^ fc fc) is the number of sections h°(X, 0^{k, k, k)) . 

By way of an example, choose the direction (1, 1, 1)Z> C if 2 (X,Z) in the Kahler 
moduli space; that is, the line bundles of the form Ojf k, k) for k G Z, k > 0. We 
list in lTable~3l the number of sections in both P 2 x P 1 x P 2 and in its restriction to the 
Schoen manifold X. Note that they grow very fast with k, and quickly grow outside 
of the range amenable to computation. However, the degree of accuracy of the metric 
on X is essentially determined by iV 2 fc k k y the number of metric parameters that we 
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fit to approximate the Calabi-Yau metric. Recall from the Hodge diamond eq. fllOOj) 
that the complex structure moduli space is 19 = /i 2,1 (X)-dimensional. However, as 
in Figure 4 , the convergence of the balanced metrics is essentially independent of the 
choice of complex structure. Hence, as an example, we choose a specific Z 3 x Z 3 
symmetric Schoen threefold (Ax = A2 = 0, A3 = 1) defined in the next section. In 



Figure 7] we plot the error measure (J(k,k,k) vs. k for this manifold and find very fast 



convergence. Note how the k = 3 data point already approaches to within 10% of 
the limit N p (= 10 6 ) > N? kkk J= 117,649), but still yields a quite small value of 

0-(3,3,3) ^4X 10" 2 . 



6 The Z 3 x Z 3 Manifold 

6.1 A Symmetric Schoen Threefold 

For special complex structures, the Schoen Calabi-Yau threefold has a free Z 3 x Z 3 
group action [4*41 151"] . which we now describe. Recall that the Schoen threefolds can 
be written as complete intersections in 

([x : x 1 : x 2 ], [t : h], [y : y x : y 2 ]) G P 2 x P 1 x P 2 , (123) 



as discussed in ISection 5l Let us start by defining the Z 3 x Z 3 group action on the 
ambient space [36], where it is generated by (u — e~ 



2tH , 
3 



and 



[xq : X\ : x 2 ] 1— > [xq : ujxi : uj 2 x-2\ 
{ [to ■ h] >-> [to : wh] (124a) 

[yo ■ yi ■ 2/2] i-> [yo ■ wyi ■ uj 2 y 2 ] 

■r„ ■ - ,■,:,,:,„ 

72 : { [t ■ h] h-> [t : *i] (no action) (124b) 

[yo ■ yi ■ 2/2] ^ [yi ■ 2/2 : yo]- 

The two generators commute up to phases on each of the two P 2 factors and, hence, 
define a Z 3 x Z 3 group action on the ambient space. Note that 72 acts non-torically, 
that is, not by a phase rotation. In order to define a Z 3 x Z 3 -symmetric Calabi-Yau 
threefold, we have to ensure that the zero locus P = = R is mapped to itself by 
the group action. For that to be the case, one must restrict the polynomials P and R 
to have a special form. It was shown in [36] that one need only constrain the cubic 
polynomials P\, P 2 , Ri, R 2 in eqns. (I102a|) and (I102bl) . Specifically, the Z 3 x Z 3 - 
symmetric Schoen Calabi-Yau threefolds are defined by the simultaneous vanishing of 
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the two polynomials 

P(x,t,y) = t P 1 (x ,x 1 ,x 2 ) + t 1 P 2 (x ,x 1 ,x 2 ) 
R(x,t,y) = t 1 R 1 (y ,y 1 ,y 2 ) + ^2(2/0, 2/i, Z/2), 



(125) 



where 



Pi (x , Xi, X 2 ) = xl + x\ + x\ + \\XqXiX 2 

P 2 (x , xt, x 2 ) = A 3 (x 2 x 2 + xjx + xlxi) 
Ri(yo, 1/1,1/2) = Vo + vt + vl + ^yovm 
R2 (2/0, 2/2) = 2/o2/i + 2/j> + vivo- 



(126) 



In the following, we will always take P, R to be of this form. Note that, up to 
coordinate changes, the polynomials depend on 3 complex parameters Ai, A 2 , and A 3 . 

One can easily check that P is completely invariant under the Z3 x Z3 group action, 
as one naively expects. However, R is not quite invariant. Rather, it transforms like 
a character of Z 3 x Z 3 . That is, 

P(jix, 7it, 7iy) =P(x, t, y) P(l2X, j 2 t, i 2 y) =P(x, t, y) (127) 

R(jix, jit, jiy) =e~^ L R(x 1 1, y) R(j 2 x, j 2 t, j 2 y) =R(x, t, y). (128) 

Nevertheless, the zero set P = = R is invariant under the group action. Moreover, 
the fixed point sets of 71 and j 2 on the ambient space P 2 x P 1 x P 2 are 

{[1:0:0], [0:1:0], [0:0:1]} X {[0:1], [1:0]} X {[1:0:0], [0:1:0], [0:0:1]}, 
{[1:1:1], [I-.w.oj 2 ], [l:w 2 :uj]j X P X {[1:1:1], [V.u:bJ 2 \, [l:u, 2 :w]}, 



respectively. For generic! 11 ! A^, the Calabi-Yau threefold X misses the Z 3 x Z 3 -fixed 
points. Therefore, the quotient 

X = X /(Z 3 x Z 3 ) = {P = = P} /(Z 3 x Z 3 ) (130) 

is a smooth Calabi-Yau threefold with fundamental group tti(X) = Z 3 x Z 3 . Its 
Hodge diamond is given by [36J 

1 


3 

h p ' q (X) =h p ' q (X/(Z 3 XZ 3 )) = 1 3 3 1. (131) 

3 

1 



11 Note, however, that Ai = A2 = A3 = is singular. A non-singular choice of complex structure 
is, for example, Ai = A2 = and A3 = 1. 
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The complex structure moduli space is h 2,l (X) = 3-dimensional and parametrized by 
Ai, A 2 , and A 3 . 



6.2 Invariant Polynomials 



As discussed in ISubsection 5.21 sections of line bundles on X are homogeneous poly- 
nomials in [xo : x\ : X2] , [to : h] and [yo : yi : J/2] , modulo the ideal (P, P) . We now 
want to consider the quotient X = X / (Z3 x Z3). Therefore, we are only interested in 
polynomials that are invariant under our group action. Let us start with the group ac- 
tion on the homogeneous coordinates (x , x±, x 2 , to, t±, y , y±, y 2 ) of P 2 x P 1 x P 2 . The 
two generators defined in eqns. ( 1124al) and ( 1124bl) can be represented by the 8x8 
matrices 
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o\ 




(0 


1 

















o\ 





to 




























1 










































1 
































1 

























1 


























UJ 











, 72 = 














1 


























1 




























1 























UJ 




























1 


\o 
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0/ 



(132) 



One can easily check that [71,72! 7^ and, in fact, the 71 and 72 actions commute up 
to multiplication by the centra] 12 ! matrix 



5 = diag(u;, uj,uj,1, 1, uj, uj, uj). 
In other words, the homogeneous coordinates 

span {x ,Xi,x 2 ,to,t 1 ,yo,yi,y2} 

= H° ( P 2 x P 1 x P 2 , 0(1, 0, 0) © 0(0, 1, 0) © 0(0, 0, r 



(133) 



(134) 



of P 2 x P 1 x P 2 carry a representation of a Heisenberg group T, which is the central 
extension 

xi><X2 

— > z 3 — ► r > z 3 x z 3 — > 0. 

Note that the map Xi x X2 is defined in terms of the two characters 

=1, 



(135) 



Xi(7i 

X2(ll) =1 



2iri 

e 3 . 



Xl(l2) 
X2(72) 



27rj 

-e 3 . 



Xi(S) =1, 
X2(5) =1 



(136) 



12 Commuting with 71 and 72. 
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of T, which will be important in the following. As discussed previously for quintics, 
ISubsection 4.21 not all line bundles are Z 3 x Z 3 -equivariant. However, computing the 
polynomials invariant under the Heisenberg group T is sufficient for our purposes. The 
T-invariants are automatically the Z 3 x Z 3 -invariant sections of Z 3 x Z 3 -equivariant 
line bundles. Their number N^ aiba ^ in each multi-degree (01,6,02) can be read off 
from the multi- variable Molien series [52] , 

P(c[x ,x 1 ,x 2 ,t ,t 1 ,y ,y 1 ,y 2 \ r ,(x,t,y)^ = £ N( aiAa2) x a H b y a2 

ai,b,a,2 

= J-£ ! 

l r l 7sr det (1-7 diag(x, x, x, t, t, y, y, y)j 

= 1 + t + t 2 + 2t 3 + 2x 3 + 2y 3 + 2x 2 y + 2xy 2 + 2t 4 + • • • . (137) 

However, to construct the Hironaka decomposition it is sufficient to determine the 
number of invariant linearly independent polynomials of total degree ai + b + a 2 . The 
corresponding Poincare series can be obtained from eq. (11371) by setting x = t = y = r, 



p(c[xo,x 1 ,X2,t ,t 1 ,yo,y 1 ,y 2 ] r ,T^ = ^iV[r fc 

k 

= ±y I 

iri ^ Hp+yi - 



r|^det(l- 7 T) 

= 1 + r + t 2 + 10r 3 + 16r 4 + 22r 5 + 85r 6 + 142r 7 + 199r 8 + 488r 9 + • • • . (138) 

Next, we need to choose 3 + 2 + 3 = 8 primary invariants. Similarly to the quintic 
case in ISubsection 4.3j we choose our primary invariants to be of the lowest possible 
degree. It is not hard to check that homogeneous polynomials 

0i = t 2 = t\ (139a) 

9 3 = XqX^ 9 4 = xl + x\ + x\ (139b) 

5 = ymV2 d & = y 3 + yl + y\ (139c) 

9 7 = x 3 xf + x 3 xl + x\x\ 8 = y 3 y 3 + y 3 y 3 + y 3 y 3 . (139d) 

can be chosen as our primary invariants. They are, in fact, the choice with the lowest 
degrees. Rewriting eq. (11381) as a fraction with the denominator corresponding to our 
choice of the primary invariants, we get 



p(c[x , xx,x 2 , t , tt, y , yi, y2\ T , r 



At 3 + 6r 4 + 6r 5 + 26r 6 + 27r 7 + 27r 8 + 46r 9 + 



(l-r)(l-r 3 ) 5 (l-r 6 )^ 

+ 42r 10 + 42T 11 + 26r 12 + 27r 13 + 27r 14 + 4r 15 + 6r 16 + 6r 17 + r 18 ) . (140) 
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deg(?7) 


# Of 7] 





1 


3 


4 


4 


6 


5 


6 


6 


26 


7 


27 


8 


27 


9 


46 


10 


42 


11 


42 


1 9 


9fi 


13 


27 


14 


27 


15 


4 


16 


6 


17 


6 


18 


1 



deg(r]) 


# of ^ 


(0,0,0) 


1 


(0,1,0) 


1 


(0,2,0) 


1 


(3,0,0) 


2 


(0,3,0) 


2 


(0,0,3) 


2 


(9 n i 


o 


(1,0,2) 


2 


(0,4,0) 


2 


(3,1,0) 


3 


(1,1,2) 


4 


(2,1,1) 


4 


(0,1,3) 


3 



Table 4 : Degrees of the 324 secondary invariants rji,... ,%24- On the left, we 
list the number of secondary invariants by total degree. On the right, 
we list some of invariants by their three individual (a±,b, 02) -degrees. 



Thus, the number of secondary invariants is 
3 5 6 2 



324 = 1 + 4 +6 + 6 + 26 + 27 + 27 + 46+ 

42 + 42 + 26 + 27 + 27 + 4 + 6 + 6 + 1. 



(141) 



Notice that the polynomials in eq. (11391) are homogeneous of multi-degree (01,6,02). 
Since the group action eq. f)132p does not mix the degrees, it follows that the secondary 
invariants will also be homogeneous polynomials. They are, moreover, separately 
homogeneous in the variables [x : X\ : x 2 ], [to : £1], and [y : yi : y 2 ]. Here, we present 
the first few secondary invariants 



>h 



(142a) 
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m 


= x 2 y yi + xiVoVz + xoVm, 




= x 2 y yi + x x y G y 2 + x yiy 2 , 


Hi 


2 2 2 

= XoVo +XiVi +x 2 y 2 , 


>k 


= x ± x 2 y + x x 2 yi + x £i2/2, 




= x 2 y + x\yi + x\y 2 , 


m 


= hyoyl + ti2/o2/2 + hyiyl, 


Vs 


= x 1 t 1 y y 1 + x t 1 y y 2 + x 2 t x y x y 2 , 




= x 2 tiyl + xohyl + xitxyl, 


Vw 


= xoX'zhyo + xoxxhy! + x x x 2 t x y 2 , 


Vu 


= xlhyo + xlhyi + xlhy 2 , 


Vl2 


= X G x\ti + x\x 2 t\ + X\x\t\, 


Vl3 


= t\ylyi + t\y\y 2 + t\y Q yl , 




= xit\yl + x 2 t\y\ + x tly 2 , 


Vl5 


= x t\y Q yi + x 2 t\y Q y 2 + x x t\y x y 2 , 



(142b) 



We list the number of secondary invariants for a given degree in ITablc 41 Thus we 
obtain the following Hironaka decomposition for the ring of T-invariant polynomials, 

324 

C[a? , xi, x 2 , t , *i, 2/o, 2/1, 2/2] r = (J) »7iC[0i, • • • , 9 8 ]. (143) 

i=l 

Finally, we need to restrict the invariant ring eq. (11431) to the complete intersection 
threefold X. In other words, one must mod out the invariant ideal 

(P, R) T = (P, R) n C[x , Xl ,x 2 , to, h, y , 2/i, 2/ 2 ] r (144) 

generated by the complete intersection equations P = = R. 

Since P is invariant, the ideal generated by P is just the invariant ring multiplied 
by P, 

324 

(P) r = 0^C[tf 1) ..,fl 8 ]. (145) 
i=i 

However, the ideal generated by R is not as simple. From eqns. (I128[) and (11361) we see 
that R transforms like the character \i- Thus the elements of the invariant ring that 
are divisible by R must also be divisible by a ^-transforming polynomial (like t\, for 
example). One can generalize the Molien formula eq. (11381) to count these "covariant" 
polynomials transforming like xi [S3], namely 

p(c[xo,Xi,x 2 ,to,ti,2/o,2/i,2/2] x? ,r N ) = -ppr , ttz . . 

V J l r l ^ det (! -7Xi(7) 2 ^) 

= r 2 + 7r 3 + 13r 4 + 22r 5 + 79r 6 + 136r 7 + 199r 8 + 478r 9 . . . (146) 
Choosing the same primary invariants as previously, eq. (11391) . one can rewrite eq. ( 11 46ft 
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as 



P( C[x , xi, x 2 , t , ti, y Q , yi, y 2 ] x \ r 



1 



(t 2 + 6r 3 + 6r 4 + 4r 5 + 27r 6 + 27r 7 + 26r 8 + 42r 9 + 



(1 -r)(l -r 3 ) 5 (l -r 6 ) 2 

+ 42r 10 + 46T 11 + 27r 12 + 27r 13 + 26r 14 + 6r 15 + 6r 16 + 4r 17 + r 20 ) . (147) 



Summing the coefficients in the numerator, we see that we again get the same number 
(= 324) of secondary x 2 -covariant generators. This is expected since we are using the 
same primary invariants. The first few secondary Xi-covariants are: 

(148a) 







77* = 




4 


= x\x 2 + xqx\ + x\x\, 


4 = 


yivl + yoyl + VoV2, 


4 


= x 2 ym + x yoV2 + xmyi, 


4 = 


x 2 yl + x y 2 + xxvl. 


4 


= x\y 2 + x\y x + x\y Q , 


4 = 


xqX-lVx + x Q x 2 y + x x x 2 y 2 , 


4 


= VoVi + v\vi + Vovlh , 


4 = 


XoXihyo + x x 2 tiy 2 + XxXitxyx 


4 


= xihvl + x 2 tiyl + x tiy 2 , 


4 = 


x 2 tiy y 2 + x tiy yi + £i*i2/i2/2, 


vi 


= x x\ti + xlxiti + x\x 2 t x , 


4 = 


x\t x y 2 + xlhyo + xltxyx, 


4i 


= x 2 t\y\ + xotjyl + x x t\y\, 


4 = 


x 2 t{yoyi + x t{yiy 2 + x 1 t 2 1 y y 2 , 



148b) 



Hence, the space of Xi-covariant polynomials, that is, transforming like Xi> is given 
by the "equivariant Hironaka decomposition" [53] (compare with eq. (11431) ) 



324 

, y 2 



C[x , x x ,x 2 , t , h, y , y u y 2 ] xl = r/f 1 C[B U ...,6 S ]. (149) 



i=i 



To summarize, even though R is not invariant, it generates an ideal which contains T- 
invariant polynomials. Using the above generalization of the Hironaka decomposition, 
a basis for these invariants is 



324 

r 



fl)=0^ 1 C[fl 1 ,. .,9s]. (150) 



i=l 



6.3 Quotient Ring 

By the results of the previous section, we know for any fixed multi-degree (a±, b, a 2 ) 
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(ai,b, a 2 ) 


N r 


N r 


(ai,6, a 2 ) 


N r 


N r 


(ai,6, a 2 ) 


N r 


N r 


(2,1,1) 


4 


4 


(2,27,1) 


56 


56 


(5,6,1) 


49 


37 


(2,2,1) 


6 


6 


(2,28,1) 


58 


58 


(5,7,1) 


56 


12 


(2,3,1) 


8 


8 


(2,29,1) 


60 


60 


(5,8,1) 


63 


47 


(2,4,1) 


10 


10 


(2,30,1) 


62 


62 


(5,9,1) 


70 


52 


(2,5,1) 


12 


12 


(2,31,1) 


64 


64 


(5,10,1) 


77 


57 


(2,6,1) 


14 


14 


(2,32,1) 


66 


66 


(5,11,1) 


84 


62 


(2,7,1) 


16 


16 


(2,33,1) 


68 


68 


(5,12,1) 


91 


67 


(2,8,1) 


18 


18 


(2,34,1) 


70 


70 


(5,1,4) 


70 


53 


(2,9,1) 


20 


20 


(3,1,3) 


23 


20 


(6,1,3) 


63 


48 


(2,10,1) 


22 


22 


(3,2,3) 


34 


29 


(6,2,3) 


94 


66 


(2,11,1) 


24 


21 


(3,3,3) 


46 


38 


(7,1,2) 


48 


38 


(2,12,1) 


26 


26 


(3,4,3) 


57 


47 


(7,2,2) 


72 


52 


(2,13,1) 


28 


28 


(3,5,3) 


68 


56 


(7,3,2) 


96 


66 


(2,14,1) 


30 


30 


(3,6,3) 


80 


65 


(8,1,1) 


30 


23 


(2,15,1) 


32 


32 


(4,1,2) 


20 


18 


(8,2,1) 


45 


31 


(2,16,1) 


34 


34 


(4,2,2) 


30 


26 


(8,3,1) 


60 


39 


(2,17,1) 


36 


36 


(4,3,2) 


40 


34 


(8,4,1) 


75 


47 


(2,18,1) 


38 


38 


(4,4,2) 


50 


42 


(8,5,1) 


90 


55 


(2,19,1) 


40 


10 


(4,5,2) 


60 


50 


(8,6,1) 


105 


63 


(2,20,1) 


42 


12 


(4,6,2) 


70 


58 


(10,1,2) 


88 


61 


(2,21,1) 


44 


44 


(4,7,2) 


80 


66 


(11,1,1) 


52 


37 


(2,22,1) 


46 


16 


(5,1,1) 


14 


12 


(11,2,1) 


78 


18 


(2,23,1) 


48 


18 


(5,2,1) 


21 


17 


(11,3,1) 


104 


59 


(2,24,1) 


50 


50 


(5,3,1) 


28 


22 


(11,4,1) 


130 


70 


(2,25,1) 


52 


52 


(5,4,1) 


35 


27 


(14,1,1) 


80 


54 


(2,26,1) 


54 


54 


(5,5,1) 


42 


32 


(14,2,1) 


120 


68 



Table 5: All homogeneous degrees leading to few (< 70 ) invariant sections 
N T = A^ lAa2) on X. For comparison, we also list the number 

N r = NF ai b a2 ) = dimC[x, t, y\Y ai & a2 ) of invariant polynomials before 
quotienting out the relations generated by the complete intersection 
equations P = = R. 
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A (finite) basis for the T-invariant polynomials 

i = c[xXy} r iaiAa2y (151) 

In particular, the polynomials are linearly independent of each other. 

Generators for the T-invariant ideal generated by the complete intersection 
eqns. f|T02a|) and ffl02bj) . 



J= (P,R) T , , , = (f>)f , , + (R) V l 

\ / (01,0,02) \ / (01,0,02) \ / (ai 



(01,6,02) \ / (01,6,02) \ / (01,6,02) 

P • C [x , t, y] T _ 3M _ ha9) , R ■ C [x , t, y] 



r ( 3bl v R.c[xXy)f bl 3) ) ■ (152) 

(01-3,6-1,02)' L ' ' WJ (»i.6-l,a2-3)/( ai)6 a2 j 



The generating polynomials of J are not automatically linearly independent. 

It remains to find a basis for the quotient 

r 

(01,6,02) 



(c[xXy]/(P,R)Y a ba = c[xXv]l lM J(P,R)l M = i/J, (153) 



corresponding to the restriction of the invariant sections on P 2 x P 1 x P 2 to the com- 
plete intersection X. This is technically more difficult than the previous quotients, 
where we were able to use Grobner bases or pick suitable primary invariants to find 
the quotient. Here, we will resort to a numerical computation of the quotient. To do 
this, note that the ideal elements J are linear combinations of invariants /. Hence, 
thinking of /, J as column vectors, there is a matrix 

M G Mat|j| X |j|(C) : MI = J. (154) 

The kernel of M is a basis for the quotient I/J. Of course, due to floating-point 
precision limits, there are generally no exact null- vectors. However, the singular value 
decomposition [51] is a well-behaved numerical algorithm to compute an orthonormal 
basis for the kernel. In ITable~5l we list the dimension 

Af aiAa2) = dime (I/J) (155) 

of the quotient space for various multi-degrees (a 1; b, a 2 ). 

6.4 Results 

We implemented Donaldson's algorithm to compute the Calabi-Yau metric on the 
threefold X = X/i/L^ x Z3). As discussed earlier, the convergence of the balanced 
metrics is essentially independent of the complex structure. Hence, we will consider 



an explicit example where Ai = A2 = 0, A 3 = 1. In Figure 8 we demonstrate that the 



numerical metric indeed approximates the Calabi-Yau metric, as it should. 
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k = k = l k = 2 k = 3 k = A 

0^(a h b, a 2 ) 



Figure 8: The error measure a^ ai ^^{X) for the metric on the Z3 x Z3- 
quotient X, computed for different Kdhler moduli but common com- 
plex structure Ai = A2 = 0, A3 = 1. Note that we chose k = 
gcd(ai, 6, 02) as the independent variable, and stopped increasing k 
as soon as N r exceeded 200. In each case we iterated the T-operator 
5 times, numerically integrating using N p = 50,000 points. Then we 
evaluated cr/ ai ^ a2 \(X) using 5,000 different test points. 
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0.1 




(iV r ) 2 = 100 (iV r ) 2 = 1000 (iV r ) 2 = 10000 
N r = H°(X,Q jt (a 1 ,b,a 2 )) T 



Figure 9: The same data as in \Figure~8 , but plotted as a function of the number 
of free parameters (-/VF ba2 ^) 2 in the ansatz for the Kahler potential. 
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In contrast to the quintic, where the single Kahler modulus is the overall volume, 
the Schoen quotient threefold X has a h 1,l (X) = 3-dimensional Kahler moduli space, 
see eq. (11311) . The Kahler moduli are determined through the three independent 
degrees (ai,b, a 2 ). Note that the integer k = gcd(ai, b,a 2 ) in |Figure~8l serves only to 
measure the refinements along a ray in the Kahler moduli space. In order to properly 
compare the metric convergence for different rays in the Kahler moduli space, we 
should consider (X^ aiba ^) 2 , which is the number of free parameters in the ansatz 
for the Kahler potential and, hence, measures the numerical complexity of the whole 
algorithm. We do this in Figure 9 and see that the accuracy is essentially determined 
by (iV£ fea2 )) 2 , and depends only slightly on the details of the Kahler moduli. 

Finally, we note again that <7( 01) & )0a )(X) is also the error measure for the metric 
pulled back to the covering space X of X^ It is useful to compare this result with the 
convergence of the Calabi-Yau metric on X obtained directly as discussed in lSection 51 
We have numerically performed this comparison and obtained results similar to those 
found in the quintic case, see Figure 4 That is, when measured by the numerical 



effort involved, the Z 3 x Z 3 symmetric method of this section is far more efficient. 
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A Primary Invariants 

In this appendix, we check that the invariants in eq. ( |88l) can be chosen to be the 
primary invariants, that is, form a "homogeneous system of parameters". In fact, the 
following criteria are equivalent, see [39] Proposition 2.3: 

• 9 2 , #3, 6*4, 6*5} are a homogeneous system of parameters (h.s.o.p.). 

• dim (C[z , z u z 2 , z 3 , z 4 ]/ (0i, 9 2 , 3 , #4, #5) ) =0 

• The only common solution to 0j = 0, % = 1, . . . , 5 is Zq = z% = z 2 = Z3 = Z4 — 0. 
Using Singular [3T], we can test the dimension criterion easily: 

SINGULAR / 
A Computer Algebra System for Polynomial Computations / version 3-0-1 

0< 
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by: G.-M. Greuel, G. Pfister, H. Schoenemann \ October 2005 

FB Mathematik der Universitaet , D-67653 Kaiserslautern \ 

> ring r=0, (z0,zl,z2,z3,z4) ,dp; 

> poly tl=z0*zl*z2*z3*z4; 

> poly t2=z0~3*zl*z4+z0*zl~3*z2+z0*z3*z4~3+zl*z2~3*z3+z2*z3~3*z4; 

> poly t3=z0~5+zl~5+z2~5+z3~5+z4~5; 

> poly t4=z0~10+zl~10+z2~10+z3~10+z4~10; 

> poly t5=z0~8*z2*z3+z0*zl*z3~8+z0*z2~8*z4+zl~8*z3*z4+zl*z2*z4~8; 

> ideal i=tl,t2,t3,t4,t5; 

> dim(std(i)) ; 




Hence, eq. ( 1551) is indeed a h.s.o.p. 

B Implementation Details 

B.l Tensors 

We use the blitz++ [55] library for all tensor computations. For example, here is the 
ordinary (serial) computation of the T-operator: 

001 template<typename Mfd> 

002 Array<C0MPLEX,2> Metric<Mf d> : : Toper at or_integr and 

003 (const Point &p) const 



003 { 




004 


using namespace blitz; 


005 




006 


Array<C0MPLEX , 1> s_val (N) , sbar_val (N) ; 


007 


s_val = X.SectionsAt(p) ; 


008 


sbar_val = conj (s_val) ; 


009 




010 


COMPLEX D = s_h_sbar(s_val, sbar_val) ; // = //s//_h"2 


011 




012 


Array<C0MPLEX,2> result (N, N) ; 


013 


firstlndex a; secondlndex b; 


014 


result = s_val(a)*sbar_val(b)/D; 


015 




016 


return(result) ; 


017 } 




018 





019 template<typename Mfd> 

020 Metric<Mfd> Metric<Mf d> : : Toperator_slow() const 

020 { 

021 Array<C0MPLEX,2> tmp(N,N); tmp = 0; 
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022 for (typename Manifold: : const_iterator p=X.begin() ; p!=X.end(); p++) 

023 tmp += Toperator_integrand(*p) * p->weight; 

024 tmp *= N / X. Volume (); 
025 

026 Metric<Manif old> result(X); result. h_ab = tmp; 

027 result . compute_hinv() ; // compute h~ab 

028 return(result) ; 

029 } 

B.2 MPI 

In order to speed up computations we use a cluster of ordinary PCs, consisting of 
11 machines connected via gigabit Ethernet. Each node has one 2.2 GHz dual-core 
Opteron processor and 2 GiB of RAM. 

The main task in computing the metrics is to compute the T-operator, see eq. ( 1271) . 
Performing the numerical integration is embarrassingly parallel, and does not even 
need any high-speed network connection. For the stochastic integration one has to 

• Compute the (weighted) integrand of the T-operator at each point, and 

• Sum the resulting matrix. 

We solve this by a bag-of-tasks, where each node adds up the contribution of a few 
points and then asks for another work set. In the end, the partial sum computed at 
each node is reduced to the master node. To effectively write distributed programs 
on the cluster, we make use of the MPI standard implemented by OpenMPI [56J. For 
example, here is the parallel implementation of the T-operator: 

001 template<typename Mfd> 

002 void Metric<Mfd>: :Toperator_slave(int) 

003 { 

004 WorkQueue<typename Mfd: : const_iterator> work (X. begin () , X.endO); 

005 Array<C0MPLEX,2> tmp(N,N) ; tmp = 0; 

006 while (work.ReceiveMoreWorkQ) { 

007 for (typename Mf d: : const_iterator 

008 p=work. begin () ; p!=work.end() ; p++) 

009 tmp += Toper at or_integr and (*p) * p->weight; 

010 } 

011 Cluster: : SendSummandArray (tmp) ; 

012 } 
013 

014 template<typename Mfd> 

015 Metric<Mfd> Metric<Mfd>: :Toperator() const 

015 { 

016 ClusterExecMethod< Metric<Mfd>, &Metric<Mf d> : :Toperator_slave > 
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017 () .Run(*this) ; 
018 

019 WorkQueue<typename Mfd: : const_iterator> work (X. begin () , X . end() , 20) ; 

020 work.FinishQ ; // main loop 
021 

022 Array<C0MPLEX,2> tmp(N,N) ; tmp = 0; 

023 Cluster: :ReceiveSumArray(tmp) ; 

024 tmp *= N / X. Volume (); 
025 

026 Metric<Manif old> result (X) ; 

027 result. h_ab = tmp; 

028 result . compute_hinv() ; 

029 return result; 

030 } 



B.3 Multivariate Polynomials 

Every section is, at the end of the day, some multivariate polynomial. For that reason, 
we implemented a C++ library for sparse multivariate polynomials. In addition to the 
usual arithmetic operations, it supports differentiation and can copy polynomials to 
remote nodes via MPI. Using this library, we can easily work with arbitrary polyno- 
mials. For example, the program to compute the metric on the Fermat quintic can, 
without change, also work with generic quintics that are a non-trivial sum over all 



126 monomials, see Figure 4 
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